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Abstract: 

Mass  transport  has  generally  been  recognized  as  a  major  limiting  factor  in  the  sensitivity  and  performance 
of  miniaturized  sensor  platforms.  To  overcome  this  limitation,  a  new  approach,  termed  molecular  dam, 
has  been  developed  to  enhance  mass  transport  for  protein  enrichment  in  nanofluidic  channels  by 
nanoscale  electrodeless  dielectrophoresis  under  physiological  buffer  conditions.  We  demonstrated  protein 
enrichment  factor  to  be  greater  than  105-fold  in  20  seconds,  which  is  orders  of  magnitude  faster  than  most 
reported  methods.  We  also  studied  the  detailed  mechanism  as  how  the  electrokinetic  forces  get  balanced 
to  achieve  the  molecular  damming  effect  and  the  scaling  relation  of  force  and  nanoconstriction  size.  Both 
the  capillary-  and  electrokinetically-driven  nanoslit  sensor  platforms  have  also  been  developed  in  parallel 
to  enhance  the  reaction  kinetics  of  protein  sensing  by  reducing  the  diffusion  length  of  reactants. 

Introduction: 

A  methodology  for  the  real-time  assessment  of  a  diverse  set  of  performance  biomarkers,  which  are 
indicators  of  fatigue,  vigilance,  and  stress,  is  currently  being  undertaken  by  the  711th  Human 
Performance  Wing  at  the  Wright  Patterson  Air  Force  Base  to  assess  the  preparedness  of  its  personnel  for 
various  missions.  However,  the  current  sensor  platform  is  unable  to  sensitively  monitor  these  biomarkers 
within  real  biological  fluid  media  such  as  saliva  and  blood,  due  to  the  presence  of  interfering  proteins  at 
far  higher  levels.  The  specific  aim  of  this  project  is  to  develop  devices  for  both  the  capillary-  (no  energy 
consumption)  and  electrokinetically-driven  nanoslit  sensors  of  relevant  neuropeptide  biomarkers,  so  that 
their  levels  can  eventually  be  monitored  in  biological  fluids.  This  project  is  to  be  performed  with  joint 
effort  in  University  of  Virginia  led  by  Prof.  Nathan  Swami  and  in  University  of  Notre-Dame  led  by  Prof. 
Hsueh-Chia  Chang.  The  collaborative  effort  with  Prof.  Swami’s  lab  is  on  the  development  of  nanofluidic 
protein  enrichment  devices  to  be  integrated  later  with  the  nanoslit  platform  for  enhancing  the  detection 
sensitivity  and  selectivity  of  biomarkers  for  human  performance  monitoring.  The  Year  1  Quarterly 
Milestones  are  listed  below  which  have  been  achieved  during  the  1st  year  of  this  project. 


Year  1  Quarterly  Milestones 

Ql 

Q2 

Q3 

Q4 

1.  Functionalization  and  fabrication  of  bio  functional  nanoslits  with  built-in 
capillary  pump 

2.  Testing  capillary-driven  nanoslit  device  with  model  protein  binding  assay 
(streptavidin-avidin)  using  fluorescence  imaging 

3.  Functionalization  and  fabrication  of  biofunctional  nanoslits  with 
electrokinetic  pump 

4.  Testing  electrokineticly-driven  nanoslit  device  with  model  protein  binding 
assay  (streptavidin-avidin)  using  fluorescence  imaging 

5.  Fabrication  of  anti-body-functionalized  nanoslits 

6.  Testing  both  electrokineticly-driven  nanoslit  device  with  antibody-antigen 
binding  assay  (IgO  anti-IgG)  using  fluorescence  imaging 

In  years  2-3,  we  envision  collaborations  with  scientists  at  the  Air  Force  Base  in  the  development  of  an 
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integrated  multi-modal  sensor  platform  for  multiple  biomarker  target  detection  for  effectively  monitoring 
the  mission  preparedness  of  its  personnel,  which  includes  initial  detection  of  Orexin  and  other  pertinent 
biomarkers  in  nanoslit  using  fluorescence  imaging,  obtain  kinetic  curves  for  titration  of  various  Orexin 
concentrations,  develop  testing  nanoslit  devices  for  spiked  samples,  and  test  nanoslit  devices  for  clinical 
samples. 

Experiments:  Description  of  the  experiment(s)/theory  and  equipment  or  analyses. 

1.  Functionalization  and  fabrication  of  biofunctional  nanoslits  with  built-in  capillary  pump 

A  single  biofunctional  nanoslit  was  fabricated  with  a  depth  of  50  nm,  followed  by  immobilizing  the 
patches  of  chemical  ligands  (biotin)  using  typical  photolithography  procedure.  The  nanoslit  was 
connected  to  periodic  structures  which  carry  the  protein  (streptavidin  conjugated  with  Alexa  Fluor  488, 
abbreviated  as  Streptavidin-488)  solution  flowing  through  the  biosensor  area.  The  immobilized  biotin 
patches  were  used  as  biosensors  to  capture  the  streptavidin-488  in  the  sample  solution. 


Figure  1.  The  schematics  of  the  biofunctional  nanoslits  with  built-in  capillary’  pump,  (a)  Streptavidin-488 
solution  was  loaded  to  the  loading  hole.  Using  the  capillaiy  pumping  structures,  streptavidin-488  were 
captured  when  aqueous  solution  flow  through  the  biotin  patterns  in  the  nanoslit.  (b,c)The  close  look  of 
the  design  of  biotin  patches  and  the  nanoslit,  (d)  The  chip  of  biofunctional  nanoslits  with  the  capillayr 
pump  structure. 


2.  Testing  capillary-driven  nanoslit  device  with  model  protein  binding  assay  (streptavidin-biotin) 
using  fluorescence  imaging. 

Streptavidin-488  was  used  as  a  model  protein  to  reveal  the  real-time  spatial  distribution  of  the  binding 
reaction.  In  figure  2,  1  nM  of  streptavidin-488  solution  was  introduced,  and  carried  into  the  nanoslit 
biosensor  without  external  pump.  Linear  flow  rate  were  estimated  to  be  25  pm/s.  Fluorescence  signal  on 
the  first  sensor  area  could  be  distinguished  after  eight  minutes,  and  reached  equilibrium  after  around  30 
minutes. 
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Figure  2.  The  real-time  binding  reaction  kinetics  of  nanoslit  biosensor  with  built-in  capillary’  pump,  (a) 
The  fluorescence  time-lapse  imaging  of  1  nM  streptavidin-488  binding  onto  the  immobilized  biotin 
patches.  Scale  bar  represents  20  um.  (b)  The  line  plot  of  the  fluorescence  intensity’  along  the  center  of  the 
nanoslits. 


3.  Functionalization  and  fabrication  of  biofunctional  nanoslits  with  electrokinetic  pump  and  testing 
electrokinetically-driven  nanoslit  device  with  model  protein  binding  assay  (streptavidin-biotin) 
using  fluorescence  imaging. 

Electrokinetic  flow  is  shown  to  be  effective  for  analyte  transport,  especially  in  nano/microfluidic  devices. 
Therefore,  we  fabricated  a  nanoslit  with  immobilized  biotin  array,  conjugated  to  microchannels  and 
external  electrodes.  For  100  nM  of  streaptavidin-488,  the  first  sensor  area  starts  to  show  significant 
fluorescence  response  after  4  minutes  and  reaches  equilibrium  after  around  6  minutes.  Using  fluorescence 
imaging,  the  spatial  distribution  of  protein  binding  events  was  recorded  in  real-time.  As  in  figure  4, 
electrokinetic  flow  increases  the  binding  events  around  -230  fold  comparing  to  simple  diffusion  in  the 
nanoslit,  and  almost  reaches  equilibrium  in  around  4  minutes.  A  notable  phenomenon  is  that  when  the 
first  sensor  area  is  reaching  its  equilibrium  fluorescence  intensity  at  around  25,000  counts,  the  second 
patch  shows  merely  around  4,000  counts,  and  the  third  to  fifth  patches  were  less  than  2,000  counts.  This 
demonstrates  the  power  of  high  capturing  efficiency  of  the  nanoslit  biosensors. 


Figure  3.  Schematics  of  the  nanoslit  biosensor  chip,  (a)  Overview  of  the  chip  setup.  Voltage  was  applied 
through  the  four  reservoirs.  Fluorescence  images  were  taken  through  the  objective  in  the  bottom,  (b) 
Three  nanoslits  were  placed  in-between  two  microchannels,  and  were  immobilized  with  biotin  patterns 
(green),  (c)  Images  taken  by  time-lapse  fluorescence  microscopy.  100  nM  of Alexa  Fluor  488-streptavidin 
was  driven  into  nanoslit  by  electroosmotic  flow  with  applied  field  of  200V/cm.  (d)  Protein  binding 
kinetics  visualized  by  fluorescence  imaging,  (insert)  100  nM  of Alexa488-streptavidin  solution  was  loaded 
and  driven  into  nanoslit  protein  sensor  patches  by  diffusion  or  electrokinetically  driven  flow.  Intensities 
in  the  center  area  along  nanoslits  (yellow  dashed  line)  was  extracted  and  plotted. 


4.  Fabrication  of  anti-body-functionalized  nanoslits  and  test  the  electrokinetically-driven  nanoslit 
device  with  antibody-antigen  binding  assay  (IgG  versus  anti-IgG)  using  fluorescence  imaging. 

For  the  antibody  immobilization  process,  the  chips  were  etched  for  microfluidic  channels  structures  and 
pre-modified  with  epoxy-silane,  which  served  as  a  linker  to  covalently  conjugate  to  the  antibodies. 
Microarray  spotter  machine  was  used  to  spatially  allocate  antibodies  onto  designed  location  on  the  chip, 
and  the  microfluidic  devices  were  encapsulated  immediately  after  the  spotting.  Chips  were  stored  in  4°C 
for  later  use.  In  figure  5,  microchannels  were  pre-wetted  with  PBS  and  later  loaded  with  PBS  containing 
fluorescence  IgG-DyFight  488  (0.15mg/mF)  and  BSA  (1.5  mg/mL)  and  driven  by  an  electric  field  of 
200V/cm.  As  shown  in  figure  5,  the  fluorescence  started  to  increase  after  30  seconds,  and  reached 
equilibrium  after  90  seconds  on  the  first  spot.  Furthermore,  all  the  three  spots  reached  their  equilibrium 
after  around  150  seconds.  Notably,  at  90  seconds,  the  second  and  third  patches  showed  significant 
darkness,  which  should  be  resulted  from  the  high  capturing  efficiency  of  the  first  spot,  which  consumes 
most  of  the  IgG-DyFight  488  flowing  through. 
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Figure  4.  Binding  kinetics  of  IgG-DyLight  488  (0.15  mg/mL)  captured 
by  the  surface  immobilized  antibodies  in  microchannel  of  1.1  \xm  depth. 
The  electric  field  was  200  V/cm,  and  the  scale  bar  represents  50  pm. 


5.  Molecular  damming  device  with  nanoconstrictions  for  ultrafast  protein  enrichment. 

The  concept  of  molecule  dam  has  been  proposed  and  device  fabricated  to  demonstrate  highly  efficient 
protein  enrichment.  Dielectric  nanoconstrictions,  down  to  30  nm  in  size,  embedded  in  nanofluidic 
device,  serve  as  field  focusing  lens  capable  of  magnifying  the  applied  field  to  105-fold  when  combined 
with  a  micro-to-nanofluidic  step  interface.  Empowered  by  this  strong  field  and  the  associated  field 
gradient  occurred  at  the  nanoconstrictions,  we  demonstrate  proteins  are  enriched  by  molecular  damming 
effect,  faster  than  the  trapping  effect,  to  greater  than  105-fold  in  20  seconds. 


Figure  5.  The  concept  of  nanoscale  molecular  trap  and  dam.  When  an  ac  electric  field  is  applied  across 
a  nanoconstriction,  composed  of  a  pair  of  sharp  insulating  tips  with  a  constriction  gap  of  10’s  nm  apart 
embedded  in  a  sealed  channel,  the  field  can  be  highly  focused  and  a  strong  field  gradient  is  generated 
for  the  operation  of  (a)  positive  dielectrophoresis  (PDEP)  or  (b)  negative  dielectrophoresis  (NDEP),  if 
the  dielectric  permittivity  of  an  analyte  (shown  in  dots)  is  larger  or  smaller  than  that  of  the  medium, 
respectively.  The  attractive  potential  in  (a)  sen’es  as  a  molecular  trap,  while  the  repulsive  potential  in 
(b)  keeps  molecules  away  from  the  constriction.  However,  if  a  dc  bias  is  applied  in  the  case  of  NDEP,  as 
shown  in  (c)  for  negatively  charged  particles,  it  tilts  the  repulsive  potential  into  a  slanted  well,  where 
the  force  balance  condition,  F ep=Feo+Fndep,  occurs  at  local  potential  minimum,  which  could  then 
cause  protein  accumulation  in  a  continuous  fashion,  effectively  working  as  a  molecular  dam. 
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Figure  6.  Demonstration  of  protein  trapping  and  damming  with  na-noscale  eDEP.  (a)  Proteins 
(Alexa-488  labeled  streptavidins,  10  pg/mL)  loaded  in  the  chip  by  capillary  force.  Protein  trapping  with 
PDEP  at  (b)  20  s  and  (c)  40  s  after  473  Vpp/cm  ac  field  applied  across  the  chip  at  10  KHz.  This  field  is 
focused  7x104  fold  (~3. 3x109  Vpp/m)  over  the  applied  field  at  the  constrictions.  All  three  constrictions 
show  trapping  of  strep-tavidins,  indicating  PDEP  is  at  work,  (d)-(f):  Protein  damming  is  demonstrated 
by  NDEP+dc:  the  ac  field  of  214  Vpp/cm  is  applied  at  1  MF4z  with  1.5  V/cm  dc  bias  (with  positive 
potential  at  the  right  side  of  the  channel).  In  (f),  the  image  was  taken  with  a  25%  transmittance  neu-tral 
density  (ND)  filter  to  keep  the  intensity  below  the  saturation  of  the  EMCCD.  The  dark  zone  at  the  center 
of  the  first  constriction  indicates  NDEP  is  at  work  where  streptavidins  are  repelled  away.  It  is  noted  that 
there  are  virtually  no  excessive  molecules  to  go  beyond  the  first  dam,  before  reaching  its  saturation,  for 
an  effective  accumulation  at  the  sec-ond  and  third  dams  around  the  constrictions  to  show  the  same  dark 
zone  as  at  the  first  constriction  with  good  contrast,  (g)  The  corresponding  intensity  plot  of  protein 
enrichment,  at  the  same  initial  concentration  of  10  pg/mL,  operated  under  the  same  experimental 
conditions  of  molecu-lar  trapping  and  damming  above.  The  plateau  in  the  damming  curve  means  the 
EMCCD  has  reached  saturation,  where  an  extended  intensity  plot  with  the  use  of  ND  filter  is  shown  in 
Figure  S2.  (h)  Protein  enrich-ment  curves  for  various  initial  concentrations  (10  ng/mL  or  189  pM,  100 
ng/mL,  and  10  pg/mL)  when  operated  under  molecular  damming  effect  with  corresponding  dc  bias  of 
1.5,  4.5,  and  1.5  V/cm.  The  concentration  rulers  of  streptavidin  are  drawn  as  horizontal  lines.  Note  that 
10  ng/mL  proteins  are  enhanced  >105  fold  in  less  than  20  seconds.  The  inset  in  (h)  shows  the  damming 
of  Alexa-488  labeled  goat  anti-human  IgG,  which  reaches  400-fold  enrichment  in  7.5  s  from  a 
concentration  of  5  pg/mL  (33  nM),  against  Alexa-488  streptavidins  (500-fold  enrichment  in  <  1  s)  under 
the  same  applied  field  conditions  as  in  (d)-(f). 


Results  and  Discussion: 

Within  nanoslits,  target  biomolecules  are  effectively  confined  to  the  surfaces,  thus  enhance  the  binding 
reactions  on  embedded  sensors.  Intuitively,  the  reduced  diffusion  time  should  shorten  the  response  time. 
However,  literatures  applying  the  protein  nanofluidic  sensors  rely  on  the  channel  electric  impedance 
measurement  display  response  times  usually  longer  than  an  hour.  To  characterize  and  optimize  the 
performance  of  nano-biosensors,  it  is  crucial  to  visualize  the  binding  process  on  sensors  in  real-time  and 
with  good  spatial  resolution,  instead  of  measuring  the  averaged  signal  changes  across  the  nanochannels. 
In  this  study,  we  investigated  the  fast  binding  kinetics  of  proteins  inside  nanofluidic  slits  due  to  the 
reduction  of  diffusion  length  dictated  by  the  slit  height.  Furthermore,  the  use  of  electrokinetic  flow  is 
shown  to  be  effective  for  analyte  transport  compared  to  pressure-driven  flow  in  nanoslits  despite  the  high 
fluid  resistance.  Our  results  present  a  platform  that  enables  visualization  of  the  reaction  kinetics  in 
nanofluidic  channels  with  potential  applications  in  kinetics  study  and  clinical  diagnostics.  In  the  2nd  year 
of  this  project,  we  aim  to  implement  the  thrombin  aptamer  sensors,  obtained  from  Air  Force  Research 
Lab,  in  the  nanoslit  to  verify  the  immobilization  protocol  and  sensing  kinetics  for  later  applications  to 
human  performance  biomarkers  detection  from  AFRL.  The  3rd  year  we  are  to  integrate  protein  enrichment 
scheme  with  molecular  dams  embedded  in  nanoslit  for  enhanced  sensing  of  these  performance 
biomarkers. 
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■  PHYSICOCHEMICAL  CLUES  TO  GENE  FUNCTION 

Messenger  RNA  (mRNA)  and  transfer  RNA  (tRNA)  are 
essential  molecules  in  protein-production  machinery:  mRNA 
serves  as  the  molecular  blueprint  that  guides  protein 
production,  and  tRNA  translates  mRNA  into  proteins  by 
joining  amino  acids  according  to  the  mRNA  sequence.  Chemist 
Bhyravabhotla  Jayaram  and  graduate  student  Garima  Khan- 
delwal  have  developed  a  computational  method  based  on 
DNA’s  physicochemical  properties  that  can  predict  whether  a 
gene  sequence  codes  an  mRNA  or  tRNA,  without  the  need  for 
structural  data  (DOI:  10.1021/ja3020956). 

The  researchers  gathered  more  than  1,500  prokaryotic 
genomes  from  public  databases,  with  more  than  2  million 
mRNA  and  50,000  tRNA  sequences.  They  used  solvation 
energy  as  a  measure  of  DNA’s  ability  to  interact  with  water  and 
calculated  for  each  genome  the  average  solvation  energy  for 
each  pair  of  neighboring  nucleotides  of  all  RNA  sequences  in 
that  genome. 

tRNA  solvation  energy  values  are  distinctly  higher  than  and 
readily  differentiated  from  those  of  mRNAs  when  comparing 
the  relative  solvation  energies  for  all  of  the  RNAs  in  each  of  the 
genomic  sequences.  This  computational  analysis  of  the 
physicochemical  characteristics  of  the  gene  sequences  reflects 
the  final  structures  for  which  they  code — tRNAs  are  stable 
structures  and  less-well  solvated  than  mRNAs.  This  solvation 
energy  analysis  offers  a  new  tactic  to  identify  tRNA-coding 
regions  in  a  genome  sequence  without  needing  tRNA’s  final 
structure,  bringing  a  new  tool  to  gene-function  studies. 
Kenneth  J.  Moore 

■  COOPERATION  YIELDS  GREAT  RETURNS 

Yonggui  Robin  Chi  and  co-workers  combine  both  N- 
heterocyclic  carbenes  (NHCs)  and  Lewis  acids,  effecting 
cooperative  catalysis  in  order  to  cleave  C— H  bonds  up  to 
three  carbons  away  from  an  aldehyde  activating  group  (DOI: 
10.1021/ja303618z).  NHC  catalysis  has  more  commonly  been 
used  to  functionalize  bonds  zero,  one,  or  two  carbons  away 
from  this  type  of  activating  group. 

Chemical  compounds  that  could  be  useful  in  synthetic, 
analytical,  or  other  practical  applications  may  have  reactive  and 
readily  functionalized  activating  groups,  as  well  as  less  reactive 
portions  of  the  molecule.  A  poorly  reactive  C— H  group,  for 
example,  may  need  to  be  functionalized  to  make  the  compound 
more  useful,  but  this  reaction  can  be  difficult  when  the  bond  is 
remote  to  the  activating  group. 

Using  an  NHC  in  combination  either  with  the  Lewis  acid 
scandium(IIl)  triflate  alone  or  with  magnesium(ll)  triflate  as 
well,  the  researchers  were  able  to  activate  a  carbon  three 
positions  away  from  the  activating  group.  After  further  reaction, 
they  constructed  6-membered  lactone  rings  that  are  important 
flavor  and  aroma  constituents  in  many  natural  products.  The 
authors  were  able  to  achieve  good  stereochemical  control,  and 
they  suggest  that  this  cooperative  catalysis  technique  might  be 
useful  for  selectively  and  asymmetrically  modifying  other 
molecules.  Christen  Brownlee 
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■  MOLECULAR  DAMS  SPEED  UP  PROTEIN 
ENRICHMENT 

Miniaturized  systems  that  rapidly  detect  minute  quantities  of 
biological  molecules  with  great  sensitivity  could  be  game- 
changers  in  medical  science.  These  micro-  and  nanoscale 
systems  have  the  potential  to  be  used  as  diagnostic  devices  to 
detect  diseases  at  preliminary  stages,  as  well  as  discovery  tools 
to  find  new  clinically  relevant  biomarkers.  But  a  major  hurdle  in 
developing  these  systems  is  picking  out  the  molecules  of 
interest  from  complex  biological  samples  and  concentrating 
them  in  one  location  for  detection. 

Now  Kuo-Tang  Liao  and  Chia-Fu  Chou  demonstrate  an 
electrokinetic  technique  that  can  quickly  concentrate  traces  of  a 
protein  at  a  single  location  (DOI:  10.1021/ja3016523).  They 
created  an  array  of  nanoscale  constrictions  that  work  like 
molecular  dams.  The  nanoconstrictions  enhance  the  local 
electric  field  by  105-fold  over  the  amount  of  applied  field. 
Proteins,  depending  on  their  size  and  dielectric  responses,  get 
trapped  at  these  nanoconstrictions. 

To  demonstrate  their  technique,  Liao  and  Chou  tested  a 
fluorescently  tagged  streptavidin  protein  molecule  and  achieved 
a  105-fold  enrichment  in  less  than  20  s.  The  investigators  say 
their  approach  is  significantly  faster  than  most  other  reported 
methods,  and  they  suggest  that  their  molecular  dams  may  be 
particularly  suitable  in  miniaturized  devices  designed  for  early 
disease  diagnostics  or  general  protein  assays. 
Rajendrani  Mukhopadhyay,  Ph.D. 

■  NEW  STUDY  REVEALS  SYNTHETIC  PRINCIPLES 
UNDERLYING  ORGANIC  SEMICONDUCTORS 

Researchers  led  by  John  Reynolds  report  synthetic  principles 
that  can  be  used  to  strategically  design  ^-conjugated  polymers 
for  applications  in  a  wide  range  of  semiconductor  technologies 
(DOI:  10.1021/ja301898h). 

Organic,  ^-conjugated  semiconductors  have  been  widely 
applied  in  recent  years  as  thin-film  transistors,  as  well  as  sensors 
and  photovoltaic  devices.  When  it  comes  to  designing  organic 
molecules  for  semiconductor  applications,  desired  character¬ 
istics  include  ease  of  processing  in  solution,  good  charge-carrier 
mobility,  a  wide  optical  absorption  profile,  and  long-term 
environmental  stability.  But  depending  on  the  desired 
application,  some  of  these  characteristics  are  more  important 
than  others.  Hence,  organic  semiconductors  designed  for  a 
specific  application  tend  to  not  be  effective  enough  for 
applications  in  other  devices,  limiting  the  material’s  broad 
usefulness. 

The  research  team  performed  a  comprehensive  study  to 
determine  the  interplay  between  the  molecular  structure  and 
the  performance  of  copolymers — made  of  alternating  units  of 
electron-rich  and  electron-deficient  repeats — in  both  thin-film 
and  photovoltaic  devices.  They  found  correlations  between  the 
charge  transport  ability  of  an  organic  semiconductor  material 
and  its  repeat  unit  structure,  molecular  weight  distribution,  and 
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ABSTRACT:  We  report  a  new  approach,  molecular  dam, 
to  enhance  mass  transport  for  protein  enrichment  in 
nanofluidic  channels  by  nanoscale  electrodeless  dielec¬ 
trophoresis  under  physiological  buffer  conditions.  Dielec¬ 
tric  nanoconstrictions  down  to  30  nm  embedded  in 
nanofluidic  devices  serve  as  field-focusing  lenses  capable  of 
magnifying  the  applied  field  to  105-fold  when  combined 
with  a  micro-  to  nanofluidic  step  interface.  With  this  strong 
field  and  the  associated  field  gradient  at  the  nano¬ 
constrictions,  proteins  are  enriched  by  the  molecular 
damming  effect  faster  than  the  trapping  effect,  to  >105-fold 
in  20  s,  orders  of  magnitude  faster  than  most  reported 
methods.  Our  study  opens  further  possibilities  of  using 
nanoscale  molecular  dams  in  miniaturized  sensing  plat¬ 
forms  for  rapid  and  sensitive  protein  analysis  and 
biomarker  discovery,  with  potential  applications  in 
precipitation  studies  and  protein  crystallization  and 
possible  extensions  to  small-molecules  enrichment  or 
screening. 


Miniaturized  biosensors  and  bioanalytical  systems  promise 
to  revolutionize  the  field  of  health  care  and  personalized 
medicine  in  light  of  sample  reduction,  speed,  and  sensitivity.1 
However,  mass  transport  has  generally  been  recognized  as  a 
major  limiting  factor  in  the  sensitivity  and  performance  of 
miniaturized  sensor  platforms,  in  that  miniaturization  leads  to 
penalties  on  passive  transport  of  biomolecules  to  the  sensor 
surface  due  to  increased  diffusion  length  from  the  bulk  liquid  or 
in  the  direction  of  the  fluidic  channels.”  Further,  for  low- 
abundance  protein  analysis,  sample  enrichment  is  often 
regarded  as  the  first  prerequisite  for  high-resolution  analysis, 
since  chemical  amplification  methods  are  not  readily  available 
for  proteins.  These  factors  pose  major  challenges  for  early  or 
acute  disease  diagnostics  and  biomarker  discovery  using  micro- 
or  nanoscale  sensor  platforms.  We  report  a  new  method  for 
rapid  enrichment  and  mass  transport  of  proteins  based  on 
electrodeless  dielectrophoresis  (eDEP)3  using  an  array  of 
insulating  nanoconstrictions  as  molecular  traps  or  dams 
(preferred  implementation),  depending  on  the  corresponding 
dielectric  response  of  the  molecules.  We  show  that  nano¬ 
constrictions,  serving  as  field  lenses,  may  enhance  local  electric 
field  to  ~105-fold  over  the  applied  ac  field  and  the  associated 
field  gradient.  As  a  result  of  this  strong  field  at  the 
nanoconstrictions,  Alexa-488-labeled  streptavidins  (52.8  kDa, 
5  nm  in  diameter)  are  enriched  by  the  molecular  damming 
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effect  to  >105-fold  in  <20  s  in  high-conductivity  physiological 
buffers,  orders  of  magnitude  faster  than  most  reported 
methods.4  Our  study  opens  further  possibilities  of  using 
nanoscale  molecular  dams  in  miniaturized  platforms  for 
sensitive  protein  analysis,  biomarker  discovery,  precipitation 
studies,  and  protein  crystallization.5 


Figure  1.  Concept  of  nanoscale  molecular  trap  and  dam.  When  an  ac 
electric  field  is  applied  across  an  insulating  nanoconstriction  with  a  gap 
of  tens  of  nm  embedded  in  a  sealed  channel,  the  field  can  be  highly 
focused,  and  a  strong  field  gradient  is  generated  for  (a)  positive 
dielectrophoresis  (PDEP)  or  (b)  negative  dielectrophoresis  (NDEP)  if 
the  dielectric  permittivity  of  an  analyte  (shown  in  dots)  is  larger  or 
smaller  than  that  of  the  medium,  respectively.  The  attractive  potential 
in  (a)  serves  as  a  molecular  trap,  while  the  repulsive  potential  in  (b) 
keeps  molecules  away  from  the  constriction.  However,  if  a  dc  bias  is 
applied  in  the  case  of  NDEP,  as  shown  in  (c)  for  negatively  charged 
particles,  it  tilts  the  repulsive  potential  into  a  slanted  well,  where  the 
force  balance  condition,  FEP  =  FE0  +  FNDEP,  occurs  at  local  potential 
minimum,  which  could  then  cause  protein  accumulation  in  a 
continuous  fashion,  effectively  working  as  a  molecular  dam. 

Among  the  efforts  for  protein  enrichment,  electrokinetic 
methods  using  ion  exclusion-enrichment  effect  by  electrical 
double-layer  (EDL)  overlapping  in  a  nanofluidic  channel,  or 
ion-selective  permeable  membranes  such  as  Nafion,  have  been 
the  most  commonly  adopted.43-1'  Other  methods  based  on 
conductivity  gradient,4’  temperature  gradient, 4k  m  and  pat¬ 
terned  membranes411  for  preconcentration  have  also  been 
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demonstrated.  However,  integrating  sensing  elements  into 
these  methods  is  not  straightforward,  since  preconcentration  is 
localized  at  the  edges  of  the  nanochannels.  On  the  other  hand, 
DEP6  and  the  recently  introduced  eDEP,3,7  or  insulator-based 
dielectrophoresis  (iDEP),8  are  known  to  be  effective  in  the 
enrichment  of  DNA,3'9  RNA,10  viruses,11  and  cells7,12  by 
exploiting  the  dielectric  response  of  the  bioparticles  and  its 
interaction  with  non-uniform  electric  fields.  In  particular,  eDEP 
and  its  integration  of  DNA  sensors  demonstrated  simultaneous 
sample  enrichment  and  sensing.1’  However,  recent  reviews  on 
eDEP  cast  doubts  as  to  whether  it  could  be  implemented  on 
chip  for  protein  enrichment  since  either  the  size  or  the 
polarizability  of  proteins  is  much  smaller  than  those  of  DNA  or 
cells,  and  hence  the  demand  for  ultrahigh  field  gradient  may  not 
be  easily  realized.14  In  the  current  study,  while  we  utilize  highly 
focused  fields  at  nanoconstrictions  to  enhance  DEP  forces,  we 
alleviate  problems  associated  with  sensor  integration  by 
applying  negative  dielectrophoresis  (NDEP)  with  dc  bias;  i.e., 
the  net  transport  from  electrophoresis  (EP)  is  opposed  by 
NDEP  and  electroosmosis  (EO)  through  the  force  balance 
condition,  Eep  =  Ee0  +  Endep,  to  enrich  proteins  away  from  the 
nanostructured  points,  a  scheme  we  call  a  molecular  dam 
(Figure  l).  This  scheme  holds  in  general  as  long  as  the  ac  field 
imposes  a  strong  NDEP  scenario  (ac  field  amplitude  dc 
bias)  at  the  constriction. 

Classical  DEP  theory  defines  the  translational  force,  i.e.,  the 
dielectrophoretic  force,  acting  on  a  polarizable  particle  in  a  non- 
uniform  field  as  Fdep  =  2sr3emTe[l,C(a))] VE2,  where  r  is  the 
radius  of  the  particle,  £m  the  absolute  permittivity  of  the 
suspending  medium,  E  the  amplitude  of  the  applied  field  (i.e., 
root-mean-squared  E  in  the  case  for  an  ac  field),  and  Ee[FC(ct))] 
the  real  part  of  the  Clausius— Mossotti  (CM)  factor, 
representing  the  frequency-dependent  dielectric  contrast 
between  the  particle  and  the  suspending  medium  in  an  external 
driving  field.4'1  It  determines  whether  the  particle  transport  is 
toward,  when  _Re[JC(co)]  >  0,  or  away  from,  when  Ee[I<C(cn)]  <  0, 
the  high  field  gradient  region  of  the  fluidic  channel, 
correspondingly  by  positive  dielectrophoresis  (PDEP)  or 
NDEP.  The  nanoscale  molecular  traps  and  dams  implemented 
in  this  study  are  depicted  in  Figure  1. 

Since  FDEP  is  proportional  to  the  size  of  the  molecules  (~r3), 
for  proteins  of  few  nanometers  in  size  (~10 — 100  kDa)  and 
small  CM  factor  due  to  low  polarizability,  it  is  challenging 
within  conventional  devices  to  enrich  proteins  by  DEP.14  To 
overcome  this  issue,  one  needs  to  create  a  highly  focused  field 
and  field  gradient  to  enhance  FDEP  by  engineering  the  VE“  (or 
E-VE)  term  described  above.  The  use  of  insulating  nano¬ 
constrictions  provides  both  an  enhanced  field  and  field  gradient 
to  compensate  the  small  size  and  low  CM  factor  of  the  proteins, 
and  to  overcome  its  large  diffusion  coefficient.  Based  on  this 
disposition,  we  have  developed  a  fabrication  process  to 
construct  eDEP  nanoconstrictions  embedded  in  nanofluidic 
channels  (nanochannels)  with  interconnections  to  microfluidic 
channels  for  sample  handling.  (Device  fabrication  is  detailed  in 
the  Supporting  Information  (SI).)  To  achieve  nanoscale  eDEP 
devices,  fused  silica  was  selected  as  the  insulating  substrate  due 
to  its  robustness  and  low  auto-fluorescence.  The  experimental 
layout  and  overall  chip  configuration  are  displayed  in  Figure 
2a, b.  Optical  and  scanning  electron  micrographs  (Figure  2c, d) 
show  that  parallel  nanochannels  containing  arrayed  nano¬ 
constriction  structures  may  be  constructed  on  one  chip. 

Due  to  the  simple  fact  that  the  displacement  current  in  a 
conducting  buffer  may  be  focused  (or  enhanced)  by  reducing  of 
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Figure  2.  (a)  Schematic  of  the  experimental  layout.  Four  Au  electrodes 
are  inserted  into  the  reservoirs,  shown  in  (b),  where  the  electric  field  is 
applied.  Experiments  are  monitored  by  an  inverted  fluorescence 
microscope  via  a  40X  objective  and  an  EMCCD  (see  Supporting 
Information),  (b)  Photograph  of  an  assembled  device  (14X14  mm2) 
made  in  fused  silica  (see  Figure  Si).  The  two  reservoirs  at  each  side  of 
the  H-shaped  microchannel,  750  pim  X  3  ftm,  are  kept  at  equal 
potential,  (c)  Optical  micrograph  showing  the  center  of  the  “H”  where 
long  nanofluidic  channels  (30  /. im  X  220  nm)  are  connected  by 
microfluidic  channels.  Three  nanoconstrictions  are  seeded  in  each  of 
the  five  nanochannels  (two  showm  here),  (d)  SEM  image  of  the  boxed 
nanoconstriction  (viewed  in  90°  turn)  in  (c),  with  30  nm  gap  size. 
Scale  bar:  (c)  30  (d)  500  nm. 


a  cross-section  of  an  insulating  fluidic  channel,  it  is 
straightforward  to  estimate  the  field-focusing  factor  in  our 
device.  If  a  microchannel  with  dimensions  -^miCroxZmicro 
(widthxheight)  is  reduced  to  a  nanochannel  of  ^nanoxZnano, 
and  the  nanochannel  is  further  reduced  to  a  nanoconstriction  of 
width  Xc,  then  from  the  continuity  equation  (current 
conservation)  one  can  easily  find  the  design  rule  for  the  overall 
field  lens  power  to  be  (Xmicro/Xnano)x(Z  micro/^nano)^  (Xnan0/ 
Xj/n,  where  n  is  the  number  of  parallel  nanochannels  (n  =  5  in 
our  device),  provided  the  conductivity  of  the  buffer  remains 
constant  over  all  fluidic  passages.  This  assumption  is  valid  as  the 
Debye  screening  length  (<1  nm)  in  the  buffers  used  in  our 
experiments  is  much  less  than  the  nanoconstriction  width  and 
nanochannel  height  of  our  device.  This  condition  is  still  far 
from  the  EDL  overlapping  scenario,  where  the  concentration 
polarization  effect  is  prounced.4a  1  From  our  current  design,  the 
cross  section  of  the  microchannel  is  750  pm  X  3  pm2,  and  that 
of  the  nanochannel  is  30  pm  X  220  nm,  with  nanoconstriction 
width  30  nm.  The  overall  field-focusing  factor  is  ~7  X  104X, 
leading  to  an  enhancement  factor  of  ~5  X  109  for  the 
dielectrophoretic  force  (~E2)  at  the  nanoconstrictions.3,7 
Empowered  by  such  a  strong  enhancement  of  the  field  and 
field  gradient,  we  were  able  to  perform  both  PDEP  and  NDEP 
in  our  device  using  Alexa-488-labeled  streptavidins  as  model 
proteins  in  physiological  buffer  conditions  (10  mM  phosphate- 
buffered  saline  with  150  mM  NaCl,  conductivity  1.6  S/m). 
Conversely,  conventional  metal  electrode-based  dielectropho¬ 
resis  (MDEP)  is  not  effective  to  overcome  the  strong  EDL 
screening  effect,  due  to  the  limitation  of  applied  field  strength, 
and  hence  could  not  easily  polarize  the  protein— counterion 
cloud  complex  in  high-salt  conditions.4a,a,h  On  the  other  hand, 
the  protein  enrichment  schemes  using  the  co-ion  depletion 
effect,  caused  by  EDL  overlapping,  suffer  from  the  relatively  low 
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potential  barrier;  hence,  high  dc  bias  cannot  be  applied  for 
rapid  protein  accumulation.4a~d,g'h 


Figure  3.  Protein  trapping  and  damming  with  nanoscale  eDEP.  (a) 
Proteins  (Alexa-488  streptavidins,  10  ^g/mL)  loaded  in  the  chip  by 
capillary  force.  Protein  trapping  with  PDEP  (b)  20  and  (c)  40  s  after 
473  Vpp/cm  ac  field  applied  across  the  chip  at  10  kHz.  This  field  is 
focused  7Xl04-fold  (~3.3Xl09  Vpp/m)  over  the  applied  field  at  the 
constrictions.  All  three  constrictions  show  trapping  of  streptavidins, 
indicating  PDEP  is  at  work,  (d— f)  Protein  damming  demonstrated  by 
NDEP+dc:  ac  field  of  214  Vpp/cm  is  applied  at  1  MHz  with  1.5  V/cm 
dc  bias  (with  positive  potential  at  the  right  side  of  the  channel).  The 
image  in  (f)  was  taken  with  a  25%  transmittance  neutral  density  (ND) 
filter  to  keep  the  intensity  below  saturation  of  the  EMCCD.  The  dark 
zone  at  the  center  of  the  first  constriction  indicates  NDEP  is  at  work 
where  streptavidins  are  repelled.  Virtually  no  excess  molecules  go 
beyond  the  first  dam  before  reaching  its  saturation  for  effective 
accumulation  at  the  second  and  third  dams  to  show  dark  zones  with 
good  contrast,  (g)  Corresponding  intensity  plot  of  protein  enrichment, 
at  the  same  initial  concentration  of  10  ^/g/mL,  under  the  same 
experimental  conditions  used  for  molecular  trapping  and  damming. 
The  plateau  in  the  damming  curve  means  the  EMCCD  has  reached 
saturation  (an  extended  intensity  plot  with  the  use  of  ND  filter  is 
shown  in  Figure  S2).  (h)  Protein  enrichment  curves  for  various  initial 
concentrations  (10  ng/mL  or  189  pM,  100  ng/mL,  and  10  /ig/mL) 
when  operated  under  the  molecular  damming  effect  with  correspond¬ 
ing  dc  bias  of  1.5,  4.5,  and  1.5  V/cm.  The  concentration  rulers  of 
streptavidin  are  drawn  as  horizontal  lines.  Note  that  10  ng/mL 
proteins  are  enhanced  >105-fold  in  <20  s.  The  inset  shows  the 
damming  of  goat  anti-human  IgG,  which  reaches  400-fold  enrichment 
in  7.5  s  from  a  concentration  of  5  //g/mL  (33  nM),  against 
streptavidins  (500-fold  enrichment  in  <1  s)  under  the  same  applied 
field  conditions  as  in  (d-f). 


Figure  3  summarizes  our  observation  of  the  molecular 
trapping  and  damming  effects  through  the  operation  of  PDEP 
and  NDEP  with  dc  bias,  respectively.  Panels  a— c  demonstrate 
that  PDEP  for  streptavidins  may  operate  in  arrayed  nano¬ 
constrictions  and  occurs  at  10  kHz.  However,  when  the 
frequency  is  increased  to  ~1  MHz,  the  DEP  undergoes  a 
crossover  response  from  PDEP  to  NDEP  with  a  changed  sign 
of  K{co),  consistent  with  other  published  results.16  Panels  d— f 
show  the  molecular  damming  effect  for  streptavidin  enrichment 
when  operated  under  NDEP  with  dc  bias.  These  observations 
verify  that  the  concepts  depicted  in  Figure  1  are  at  work. 

To  compare  the  effectiveness  of  protein  enrichment  by 
PDEP  (trapping)  and  NDEP  with  dc  bias  (damming),  we 
loaded  the  chips  with  the  same  initial  concentration  of  10  /,<g/ 


mL  Alexa-488  streptavidins.  Figure  3g  indicates  that  the 
damming  effect  is  much  more  efficient  in  protein  enrichment, 
as  103-fold  concentration  enhancement  may  be  achieved  2—3  s 
after  the  field  is  turned  on,  much  faster  than  the  trapping  effect 
(see  SI  Figure  S2  and  Movies  Ml  and  M2).  In  fact,  this  result 
suggests  that,  for  practical  applications,  nanoscale  eDEP  as  a 
molecular  dam  (NDEP+dc)  is  particularly  advantageous  for 
three  reasons:  it  is  much  more  effective  in  protein  enrichment 
than  PDEP;  concerns  of  potential  Joule  heating  effect  are 
alleviated  by  displacement  of  the  molecular  dam  away  from  the 
hot  spot,  i.e.,  the  geometrical  center  of  the  nanoconstriction; 
and  the  sensing  element  can  be  placed  micrometers  away  from 
the  nanoconstriction,  a  task  easily  achievable  by  conventional 
photolithography.  This  circumvents  the  great  technical 
challenge  of  integrating  biosensors  at  the  center  or  in  the 
immediate  proximity  of  the  nanoconstrictions  when  operating 
in  PDEP  scheme. 

To  further  characterize  the  enrichment  factors  for  different 
initial  protein  concentrations  under  molecular  damming 
conditions,  streptavidins  of  10  ng/mL,  100  ng/mL,  and  10 
fig/mL  were  loaded  into  the  chips.  Figure  3h  shows  the 
enrichment  curves  at  various  starting  concentrations  of 
streptavidin,  where  the  concentration  rulers  of  1  and  5  mg/ 
mL  are  drawn  as  horizontal  lines  under  the  same  buffer 
conditions  as  those  used  in  experiments  to  serve  as  references 
for  calculating  the  enrichment  factor  and  the  time  needed  to 
reach  the  concentration  rulers.  To  ensure  the  molecular 
damming  strategy  developed  here  is  applicable  to  other 
proteins,  we  also  tested  Alexa-488-labeled  goat  anti-human 
IgG  (~150  kDa,  see  SI)  in  comparison  to  Alexa-488 
streptavidins  (inset  in  Figure  3h).  The  results  indicate  the 
effect  of  the  molecular  size  and  the  polarizability  for  different 
proteins.  Note  that  all  the  data  presented  in  Figure  3g,h  and 
Figure  S2  are  plotted  from  the  region  of  highest  intensity  of 
fluorescent  signals  (with  an  area  of  9.6  //m2  or  24  pixels)  where 
proteins  are  mostly  enriched,  after  subtracting  the  background 
from  dark  counts  and  auto-fluorescence  from  the  substrate. 
Based  on  these  results,  a  protein  enrichment  factor  of  >105-fold 
may  be  achieved  in  just  seconds  when  operating  in  the 
damming  scheme,  orders  of  magnitude  faster  than  most 
previous  studies.4a,c,h  We  attribute  the  fast  transport  of 
molecules  to  the  highly  constricted  field  at  both  the  micro- 
to  nanochannel  junction  and  the  nanoconstriction.  This  is 
indeed  a  unique  feature  of  our  device.  Streptavidin  velocity  is 
enhanced  from  ~1.5  /.im/s  in  microchannel  (1.5  V/cm  applied 
dc  bias  and  bulk  mobility  0.8  ±  0.9  ^fm-cm/V-s)4'  to  ~100  [A m/ 
s  by  a  70-fold  field  enhancement  when  entering  into  the 
nanochannel,  and  further  to  ~10— 15  cm/s  by  another  1000- 
fold  field  enhancement  from  the  junction  to  the  nano¬ 
constriction  (~45  fim  distance).  To  estimate  the  effective 
potential  energy  Umin  involved  in  the  damming  process  using 
Boltzmann  distribution,  105-fold  concentration  enrichment 
corresponds  to  Umin  «  — 12  kBT,  where  kBT  is  the  thermal 
energy.  Since  in  the  experiments  the  system  did  not  reach 
equilibrium  as  proteins  kept  accumulating  during  the  process, 
this  estimation  would  suggest  a  lower  bound  of  the  effective 
potential  well. 

To  demonstrate  that  the  balance  of  FEO  and  f^EP  alone  cannot 
achieve  the  significant  protein  enrichment  observed  in  our 
device,  we  studied  the  case  of  a  pure  dc  field  (up  to  4.5  V/ cm) 
but  did  not  observe  any  discernible  protein  enrichment  (see  SI, 
Movie  M3).  In  fact,  it  suggests  FEP  is  much  higher  than  FE0  in 
our  device.  To  further  demonstrate  the  essential  role  of  FNDEP 
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in  molecular  damming,  we  turned  off  the  ac  field  but  not  the  dc 
bias  and  observed  the  highly  enriched  proteins  quickly  diffusing 
away  due  to  the  concentration  gradient  established  by  the 
molecular  dam  (see  SI,  Movie  M4).  Hence,  TNDEP  is  essential 
for  the  protein  damming  effect  in  our  design. 

The  potential  Joule  heating  effect  due  to  the  highly  focused 
field  (total  current  25  \Jt A ,  or  5  //A/nanochannel)  at  the 
nanoconstrictions  may  be  alleviated  by  the  small  sample 
volume  (~1  pL/nanochannel),  with  a  220  nm  liquid  layer  in 
the  nanochannels,  used  in  our  devices.  Hence  heat  dissipation 
through  the  substrate,  as  a  bulk  thermal  bath,  is  very  effective 
based  on  our  previous  study  by  finite-element  multiphysics 
simulation.13'17  Within  our  experimental  conditions,  we  assume 
the  proteins  at  the  trap  are  not  denatured,  as  the  trapping 
events  are  reversible,  and  there  is  no  denaturation-associated 
aggregation  observed  in  our  experiments.  The  issues  discussed 
here  are  further  assuaged  if  one  operates  the  device  as 
molecular  dams  where  enriched  proteins  are  away  from  the 
highly  focused  field  constrictions. 

Although  protein  trapping  has  been  demonstrated  using  a 
100-nm  nanopipet  with  quasi-dc  driving  field,18  and  insulator 
post-array  with  dc  field  (iDEP  with  low-conductivity  buffers  to 
avoid  Joule  heating),19  it  remains  challenging  to  detect  ultra- 
low  protein  targets  due  to  the  limited  protein  enrichment  factor 
achieved  (up  to  1000-fold)  and  to  integrate  it  within  a  fluidic 
chip  platform  where  multiplexing  and  parallel  analysis  are 
desirable  (in  the  nanopipet  approach).  However,  the  concept  of 
the  molecular  damming  effect  introduced  in  this  study  is 
exceedingly  compatible  with  multiplexing,  parallel  analysis,  and 
high-conductivity  buffers,  and  thus  suitable  for  integrating 
biosensors. 

In  summary,  we  present  a  nanoscale  active  molecular 
transport  scheme  for  ultrafast  protein  enrichment  by 
constructing  nanoscale  molecular  traps  and  dams  using 
electrodeless  dielectrophoresis  generated  by  insulating  nano- 
constricted  structures.  A  protein  enrichment  factor  >105  has 
been  achieved  in  <20  s,  orders  of  magnitude  faster  than  most  of 
the  reported  methods.  Multichannel  layout  for  parallel 
operation  has  also  been  demonstrated.  In  this  scenario, 
miniaturization  alleviates  rather  than  accentuates  the  transport 
limitations,  so  that  any  sensor  applications  can  capitalize  on  the 
ultrafast  sample  enrichment  schemes  introduced  here.  Though 
our  device  could  also  be  applied  to  DNA  and  RNA 
analysis,3’10'13  it  may  find  applications  in  general  protein  assays, 
protein  crystallization,  rare  biomarker  discovery  (e.g.,  coupled 
with  mass  spectroscopy),  and  early  disease  diagnostics  in  lab- 
on-a-chip  systems,  with  potential  extensions  to  enrichment  or 
screening  of  small  molecules  (e.g.,  peptides  or  carbohydrates). 
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Research  Article 

Nano-constriction  device  for  rapid  protein 
preconcentration  in  physiological  media 
through  a  balance  of  electrokinetic  forces 

We  describe  a  methodology  to  steeply  enhance  streptavidin  protein  preconcentration 
within  physiological  media  over  that  achieved  by  negative  dielectrophoresis  (NDEP) 
through  utilizing  a  DC  offset  to  the  AC  field  at  nanoscale  constriction  gap  devices.  Within 
devices  containing  approximately  50-nm  constriction  gaps,  we  find  that  the  addition  of 
a  critical  DC  field  offset  (1.5  V/cm)  to  the  NDEP  condition  (~200  Vpp/cm  at  1  MHz) 
results  in  an  exponentially  enhanced  extent  of  protein  depletion  across  the  device  to  cause 
a  rapid  and  steeply  rising  degree  of  protein  preconcentration.  Under  these  conditions, 
an  elliptical-shaped  protein  depletion  zone  that  is  extended  along  the  device  centerline 
axis  forms  instantaneously  around  the  constrictions  to  result  in  protein  preconcentration 
along  the  constriction  sidewall  direction.  Through  a  potential  energy  diagram  to  describe 
the  electrokinetic  force  balance  across  the  device,  we  find  that  the  potential  energy  barrier 
due  to  NDEP  is  gradually  tilted  upon  addition  of  DC  fields,  to  cause  successively  steeper 
potential  wells  along  the  sidewall  direction  for  devices  containing  smaller  constriction 
gaps.  Hence,  for  approximately  50-nm  constriction  gaps  at  a  critical  DC  field,  the  ensuing 
narrow  and  deep  potential  energy  wells  enable  steep  protein  preconcentration,  due  to 
depletion  over  an  exponentially  enhanced  extent  across  the  device. 


Keywords: 
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1  Introduction 

The  challenge  of  sensing  rare  numbers  of  biomarker  pro¬ 
teins  against  a  background  of  high  concentration  of  other 
matrix  proteins  within  physiologically  relevant  media  re¬ 
quires  methodologies  for  selective  preconcentration  of  the 
biomarker  in  the  proximity  of  the  sensor  [1],  Conventional 
chemical  methods  based  on  antibody  depletion  are  unable 
to  achieve  the  necessary  degree  of  preconcentration,  since 
the  biomarkers  are  present  at  106-1012-fold  lower  levels 
than  the  background  proteins  in  blood  [2],  Hence,  there  is 
great  interest  in  applying  electrokinetic  methods,  especially 
within  nanofluidic  devices  [3],  where  the  required  degrees  of 
preconcentration  may  be  achieved  due  to  the  large  volume 
reduction.  While  prior  approaches  have  investigated  pre¬ 
concentration  using  ion  exclusion-enrichment  effects  due 
to  electrical  double  layer  overlap  at  the  micro-to-nanofluidic 
interface  [4-6],  the  resulting  protein  preconcentration  at 
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103-106-fold  levels  requires  several  tens  of  minutes  to  an 
hour,  and  this  method  cannot  effectively  separate  target 
biomarker  proteins  against  a  background  of  similarly  charged 
and  sized  proteins  in  samples  of  physiological  fluids.  Di¬ 
electrophoresis  (DEP)  enables  highly  selective  trapping  of 
bioparticles  based  on  the  characteristic  frequency  response 
of  the  dielectric  permittivity  of  the  bioparticle  versus  that 
of  the  medium  [7,  8],  and  it  has  been  extensively  applied 
toward  sorting  of  somewhat  similar  sized  biological  cells 
with  differing  dielectric  frequency  response  [9].  However,  its 
application  to  smaller  sized  biomarlcers,  such  as  nanoscale 
proteins  and  fragments  of  ss-DNA,  requires  methods  to  en¬ 
hance  the  local  field  to  offset  the  steep  fall  in  DEP  trap¬ 
ping  forces  with  particle  size  [10, 11].  Recently,  through  the 
application  of  devices  with  500-fold  dielectric  constrictions 
(500  to  1  p,m)  to  enhance  the  field,  we  have  demonstrated 
the  trapping  of  150-base  ss-DNA  in  media  of  high  con¬ 
ductivity  to  enable  target  DNA  preconcentration  and  sens¬ 
ing  through  hybridization  with  immobilized  capture  probes 
[12-14].  However,  given  the  smaller  sizes  of  proteins,  the  di¬ 
rect  trapping  of  proteins  by  DEP  requires  much  higher  fields. 
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Hence,  higher  voltages  are  needed  for  trapping  within  an 
electrode  device  format  [15,16],  which  can  affect  protein  func¬ 
tionality,  and  within  electrode-less  configurations,  it  requires 
nano-sized  structures,  such  as  proof-of-concept  devices  uti¬ 
lizing  nano-pipettes  [17]  or  nanotubes  tips  [18],  which  are 
not  scalable.  Other  trapping  methods  within  droplets  [19], 
pose  limitations  on  the  media  for  protein  preconcentration. 
Alternatively,  indirect  means  can  be  applied  to  manipulate 
proteins  through  streaming  DEP  [20,21],  to  influence  the  net 
flow  profiles  in  a  microfluidic  device,  but  the  degree  of  control 
of  protein  trapping  to  enable  preconcentration  and  separation 
is  limited  by  the  small  magnitude  of  dielectrophoretic  forces 
over  other  forces.  Hence,  optimizing  the  dielectrophoretic 
trapping  of  proteins  within  a  scalable  nanofluidic  device  is  of 
interest,  to  assess  the  influence  of  device  structure  (surface 
charge  at  constriction  tip  and  local  electric  field),  electrical 
conditions  (AC  bias  frequency  and  DC  bias),  and  fluidic  con¬ 
ditions  (media  conductivity  and  pH)  on  protein  preconcen¬ 
tration  and  separation. 

The  trapping  of  proteins  at  room  temperature  condi¬ 
tions,  within  physiological  media  of  neutral  pH  and  high 
conductivity  is  necessary  to  maintain  the  conformation  and 
biofunctionality  of  the  protein.  However,  DEP  trapping  of 
proteins  within  media  ofhigh  conductivity  poses  distinct  chal¬ 
lenges.  First,  the  frequency  response  of  the  protein  species 
in  conducting  media  is  difficult  to  interpret  given  the  com¬ 
plex  interplay  of  surface  charges,  mobile  charges,  and  water 
molecules  under  the  electric  field,  especially  since  counter¬ 
ion  fluctuations  influence  the  relaxation  time  for  polariza¬ 
tion  of  biomolecules  [22,  23].  Second,  the  enhanced  current 
flow  due  to  localized  fields  within  media  of  high  conduc¬ 
tivity  (>1  S/m)  causes  Joule  heating  [24,  25],  which  can 
denature  the  proteins.  Finally,  while  this  problem  can  be 
overcome  by  applying  negative  DEP  to  trap  proteins  away 
from  the  high-field  region,  the  resulting  temperature  gra¬ 
dients  due  to  Joule  heating  can  cause  bulk  electrothermal 
flow  that  can  lower  DEP  preconcentration  [13].  Hence,  in 
order  to  enhance  the  degree  of  preconcentration,  there  is  a 
need  of  methods  to  localize  protein  trapping  under  negative 
DEP  (NDEP).  Herein,  through  the  application  of  a  scalable 
nanofluidic  device  platform,  we  investigate  the  influence  of 
DC  bias  coupled  to  the  AC  field  on  devices  of  successively 
smaller  constriction  gaps,  to  enable  narrowing  of  the  protein 
preconcentration  band  through  a  balance  of  electrophoresis, 
electro-osmosis,  and  NDEP  forces.  We  find  that  steeply  ris¬ 
ing  levels  of  protein  preconcentration  can  be  obtained  only 
under  conditions  that  enable  a  depletion  zone  of  enhanced 
extent,  which  occurs  almost  instantaneously  in  proximity  of 
approximately  50  nm  constriction  gaps,  at  frequencies  char¬ 
acteristic  of  NDEP  (1  MHz),  and  in  presence  of  a  critical 
DC  bias  (~1.5  V/cm).  Using  a  potential  energy  diagram  to 
describe  the  electrokinetic  force  balance  across  the  device, 
we  find  that  application  of  approximately  50  nm  constric¬ 
tions  to  enable  a  significant  energy  barrier  to  protein  trans¬ 
port,  in  conjunction  to  a  critical  DC  field  to  cause  significant 
tilting  of  the  energy  profile  results  in  the  formation  of  nar¬ 
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row  and  deep  potential  wells  for  rapid  protein  preconcentra¬ 
tion,  due  to  protein  depletion  over  an  exponentially  enhanced 
extent  versus  devices  under  subcritical  DC  or  with  larger 
constrictions. 


2  Materials  and  methods 

2.1  Nanoconstriction  device 

The  fabrication  of  nanofluidic  devices  with  nanoconstrictions 
for  ultrafast  protein  preconcentration  is  based  on  our  previ¬ 
ously  developed  insulator  constriction-based  dielecrophore- 
sis  device  (see  Supporting  Information)  [26].  A  microfluidic 
channel  was  constructed  by  photolithography  on  a  fused- 
silica  substrate,  of  3  p.m  depth  and  750  p,m  width,  with  outer 
constrictions  to  form  five  “nanochannels”  of  200  nm  depth 
and  30  p,m  width.  Following  this,  electron  beam  lithogra¬ 
phy  was  applied  to  nanofabricate  multiple  inner  constrictions 
within  the  nanochannel  of  varying  sizes  (15  nm  to  150  nm), 
as  confirmed  by  scanning  electron  microscopy.  The  whole 
device,  with  integrated  micro  (3  p,m  depth)  and  nanofluidic 
(200  nm  depth)  channels  was  completed  at  room  temper¬ 
ature  using  a  low-pressure  sealing  process  [27].  The  final 
nanofluidic  device  is  shown  in  Fig.  SI  of  the  Supporting 
Information. 


2.2  Protein  preconcentration  and  imaging 

This  device  was  subsequently  applied  toward  rapid  protein 
trapping  under  positive  DEP  (PDEP)  at  approximately  200 
Vpp  across  approximately  1  cm  at  100  kHz,  and  negative  DEP 
(NDEP)  at  approximately  200  Vpp  across  approximately  1  cm 
at  1  MHz,  using  an  Agilent  function  generator  (33220A)  and 
an  FLC  voltage  amplifier  (A400D).  To  enhance  preconcen¬ 
tration  under  NDEP,  we  investigated  additional  DC  offsets 
of  approximately  0.3  V/cm  and  approximately  1.5  V/cm  field 
strength.  The  preconcentration  of  fluorescently  labeled  strep- 
tavidin  protein  samples  (52.8  kDa,  Molecular  Probes,  Engene, 
OR,  USA),  at  p,g-ng/mL  concentration  levels,  within  an  elec¬ 
trolyte  of  10  mM  PBS  containing  150  mM  NaCl  and  2  mM 
NaN3,  at  pH  7.2,  was  imaged  using  an  inverted  Zeiss  Ob¬ 
server  microscope  with  a  Hammatsu  EM-CCD  to  enhance 
sensitivity. 


2.3  Modeling  methods 

The  electric  potential  within  the  nanofluidic  device  can  be 
mapped  through  solving  the  Laplace  equation  to  determine 
the  equipotential  lines: 

V  •  ((cr  +  j  coe)  V<j>)  =  0  (1) 
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The  electric  field  distribution  can  then  be  obtained  using 
the  relation: 

E  =  -V<|>  (2) 

Based  on  the  electric  field  distribution  in  the  device, 
its  gradient  can  be  calculated,  which  in  turn  can  be  used 
to  determine  the  dielectrophoretic  force  profile  across  the 
device.  At  steady  state,  the  average  dielectrophoretic  force, 
Fdep,  on  a  homogeneous  spherical  particle  with  electri¬ 
cal  permittivity,  ep,  conductivity,  ap,  and  radius,  a,  sus¬ 
pended  in  a  fluid  with  electrical  permittivity,  em,  and  con¬ 
ductivity,  o-m,  in  a  nonuniform  AC  electric  field,  E,  is 
given  by  the  following  in  3D  with  a,  representing  the 
polarizability: 

- —  |V£2  =  a£.V£  (3) 

+  2e*M  J  ' 

- v - ' 

kcm 

Here,  a  =  2Ttsma3  (  .e'"  )  and  KCm  is  the  Clausius- 

\  £p-r£Em  / 

Mossotti  factor,  which  depends  on  the  frequency  (o>)  of  the 
applied  field  as  given  by:  e*  =  e  +  a/joo.  Since  the  DEP 
force  is  experienced  by  the  polarized  particles,  only  in  nonuni¬ 
form  electric  fields,  it  is  proportional  to  the  gradient  (or  spa¬ 
tial  nonuniformity)  of  the  square  of  the  electric  field.  The 
electrophoresis  force  (FEP)  can  be  calculated  from  the  elec¬ 
trophoretic  mobility  (|xEp),  using  friction  coefficient  (y  = 
6rrqa;  for  protein  particles  of  size,  a,  and  medium  viscosity: 
iq)  as: 

Fep  =  yp.EPE  (4) 

Similarly,  the  electro-osmosis  force  (FE0)  can  be  calcu¬ 
lated  from  the  electro-osmotic  mobility  (|_lep)  within  the  de¬ 
vice  and  friction  coefficient  (y): 

Fbo  =  TM-eoE  (5) 

As  in  Fig.  1,  under  NDEP  conditions,  FNDEP  and  FE0 
oppose  FEp  on  the  side  of  the  constriction  facing  the  inlet,  as 
given  by  net  force  (Fnet): 

Fnet  =  Fep  —  Feo  —  Fndep  (6) 

On  the  side  of  the  constriction  facing  the  outlet,  where  no 
protein  trapping  occurs,  FNDEP  and  FEP  oppose  FE0,  thereby 
reversing  the  sign  on  FNDEP  in  Eq.  (6).  Since  all  steady-state 
electrical  forces  are  conservative  vector  fields,  we  can  calculate 
scalar  potential  energy  fields,  to  describe  the  net  electrokinetic 
force  balance  in  terms  of  barriers  and  wells  that  affect  particle 
trapping  within  the  device.  The  potential  field  is  given  by  the 
integral  over  the  volume  (dV)  of  the  net  electrokinetic  force 
(Fnet): 

U=j)Fnet.dV  (7) 


Based  on  Eqs.  (3-6),  this  can  be  written  as: 

U=  (f  aEVEdV+  (f  y(|xep  +  |xe0)  •  EdV  (8) 

J  v  Jv 

This  can  be  integrated  to  the  following,  based  on 
the  integral  forms  of  Eq.  (2)  for  the  DC  voltage, 
$DC: 

aE 2 

u  =  — - y(tcp  +  |XE0)  •  ^DC  (9) 


3  Results 

3.1  Device  methodology  for  protein 
preconcentration 

The  methodology  for  the  application  of  nanoscale  constric¬ 
tions  to  enable  protein  preconcentration  through  a  balance 
of  electrokinetic  forces  is  schematically  described  in  Fig. 
1.  Under  NDEP  conditions  (~200  Vpp/cm  at  1  MHz),  pro¬ 
tein  molecules  are  pushed  away  from  the  constriction  edge 
by  FNDEP,  as  shown  in  Fig.  1A.  Within  physiological  media 
(~150  mM  salt  at  pH  7.2)  of  high  conductivity  (am  ~1.6 
S/m),  the  enhanced  current  and  field  cause  a  substantial 
level  of  Joule  heating  (see  Supporting  Information  Fig.  S3), 
with  the  resulting  temperature  gradients  causing  a  net  elec¬ 
trothermal  flow.  Based  on  prior  work  [28],  the  direction  of 
electrothermal  flow  (FETh)  is  given  by  -VT-E.  Within  con¬ 
striction  DEP  devices,  since  the  temperature  gradient  (VT) 
and  the  electric  field  (E)  vectors  are  pointed  toward  the  con¬ 
striction  edge  [13],  FEth  is  pointed  away  from  the  constriction 
edge,  similar  to  the  direction  of  FNDEP.  However,  FNDEP  is 
highly  localized  at  the  constriction  due  to  its  dependence  on 
VE2,  whereas  FETH  is  more  long  ranged  [24,  25],  since  it  de¬ 
pends  on  E2  and  acts  indirectly  upon  the  particles  through 
a  drag  force.  Hence,  under  NDEP  conditions  within  media 
of  high  conductivity,  the  protein  trapping  is  spread  over  a 
diffuse  region,  thereby  requiring  a  DC  field  to  enable  the 
observation  of  protein  preconcentration.  Upon  the  applica¬ 
tion  of  a  small  DC  field  in  addition  to  the  higher  AC  bias  at 
1  MHz,  we  enable  the  action  of  electrophoretic  forces  (FEP)  on 
the  streptavidin  protein  particles,  as  well  as  electro-osmotic 
drag  forces  (FEo)  from  the  fluid,  in  conjunction  to  the  NDEP 
forces.  As  shown  in  Fig.  IB,  under  electrophoresis,  the  nega¬ 
tively  charged  protein  molecules  (negative  2.3e,  as  per  [3]  for 
an  elementary  charge  of  “e”)  are  drawn  away  from  the  termi¬ 
nal  at  negative  potential  (henceforth  called  "inlet”) ,  toward  the 
nanoconstriction.  On  the  other  hand,  under  electro-osmosis 
and  negative  DEP  forces,  the  molecules  are  pushed  away  from 
the  nanoconstriction  and  back  toward  the  inlet.  In  the  vicin¬ 
ity  of  the  1000-5000-fold  constrictions  used  within  this  work, 
Fndep  is  substantially  larger  than  FEP  and  FE0.  Since  protein 
preconcentration  and  depletion  only  occur  on  the  side  of  the 
constriction  facing  the  inlet,  we  assume  that  |xEP  exceeds  p,Eo 
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Figure  1.  Schematics  of  protein  pre¬ 
concentration  through  an  electrokinetic 
force  balance:  (A)  Trapping  under  NDEP 
(NDEP)  using  AC  bias  only;  (B)  Apply¬ 
ing  an  additional  DC  offset  over  the  AC 
bias  can  localize  protein  preconcentra¬ 
tion  between  the  inlet  and  the  constric¬ 
tion,  as  shown  in  (C)  through  a  balance  of 
electro-osmosis  (Feo)  and  NDEP  (Fndep) 
versus  electrophoresis  (Fep)  to  cause  for¬ 
mation  of  the  "depletion  zone"  of  arc  ra¬ 
dius,  r,  near  the  constriction  (grey),  and 
the  "preconcentration  zone"  of  width  w. 
The  effect  of  constriction  size  (d)  and  ap¬ 
plied  field  conditions  (DC  versus  AC  field) 
on  protein  preconcentration  is  studied  by 
investigating  the  force  balance  along  the 
centerline  (V)  and  sidewall  vectors  ( Y ). 


at  all  points  across  the  device,  so  that  in  the  proximity  of  the 
constriction  where  Fndep  is  substantial,  the  sum  of  electro¬ 
osmosis  and  NDEP  forces  can  just  equal  the  opposing  elec¬ 
trophoresis  force,  as  given  by:  FEP  =  FE0  +  FNDEP.  This  would 
explain  protein  preconcentration  onwards  from  this  force  bal¬ 
ance  point,  with  protein  depletion  occurring  in  the  vicinity  of 
the  constriction,  as  in  Fig.  1C.  This  assumption  of  |xEP  >  p,E0 
also  explains  the  absence  of  protein  depletion  and  preconcen¬ 
tration  on  the  side  of  the  constriction  facing  the  outlet  where 
Fndep  opposes  FE0.  Fience,  the  force  balance  condition  of: 
Fep  +  Fndep  =  FEo  will  not  be  reached,  given  that  FEP  >  FE0, 
and  the  addition  of  FNdep  to  FEP  would  continue  to  be  greater 
than  FEo-  The  high-salt  levels  within  the  medium  (~150  mM 
salt)  serve  to  enhance  |xEP  over  p,E0,  since  the  thinner  electri¬ 
cal  double  layer  around  the  channel  wall  reduces  (xE0  and  the 
smaller  counter-ion  cloud  around  protein  molecules  reduces 
retardation  forces  to  electrophoresis  arising  from  hydrody¬ 
namic  drag  and  ion  flow  within  the  counter-ion  cloud  [29]. 
We  have  assumed  p,EP  of  approximately  1  p,m.cm.V_1s_1  [3]; 
and  approximately  10-fold  lower  p,Eo,  with  the  DEP  mobility 
making  up  the  difference  in  proximity  of  the  constriction,  as 
computed  from  an  NDEP  force  using  a  Clausius-Mossotti 
factor  (KCm)  of  —0.5.  These  depletion  and  preconcentration 
zones  occur  along  an  arc,  as  in  Fig.  1C.  In  this  manner,  the 
trapping  region  can  be  narrowed  to  enhance  the  degree  of  pro¬ 
tein  preconcentration  in  comparison  to  what  can  be  achieved 
using  NDEP  only.  The  degree  of  protein  preconcentration  can 
be  further  enhanced  by  the  development  of  conditions  to  en¬ 
hance  the  extent  of  the  depletion  zone  (r)  or  reduce  the  width 
of  the  preconcentration  zone  (w).  In  the  following  sections,  we 
investigate  the  critical  values  of  the  DC  field,  AC  frequency, 
and  constriction  size  (d)  required  to  cause  rapid  enhancement 
of  protein  preconcentration  through  exploring  the  profile  of 
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the  electrokinetic  force  balance  (Fnet,  in  Eq.  (6)  across  the 
device.  Specifically,  we  compare  the  force  balance  along  di¬ 
rections  of  the  centerline  (vector  Y  along  device  length)  and 
sidewall  (vector  Y'  along  constriction  edge),  as  in  Fig.  1C, 
to  track  formation  of  the  depletion  and  preconcentration 
zones. 

To  optimize  the  electrokinetic  force  balance  for  enhanced 
protein  preconcentration,  the  time  evolution  of  protein  pre¬ 
concentration  was  investigated  at  devices  of  varying  nanocon¬ 
striction  gap  size  (30  p,m  channel  constricted  to  d  ~50  nm, 
~100  nm  or  ~140  nm)  and  under  varying  DC  field  offsets 
(~0.3  V/cm  or  ~1.5  V/cm)  over  the  AC  field  conditions 
for  NDEP.  Figure  2  shows  time-lapse  fluorescence  images 
for  protein  trapping  in  the  region  close  to  the  device  con¬ 
striction  under  NDEP  conditions  (~200  Vpp/cm  AC  field  at 
1  MFIz)  with  an  added  DC  field  (~1.5  V/cm)  for  devices  with 
constrictions  of  approximately  50  nm  (Fig.  2A-C)  versus  ap¬ 
proximately  140  nm  (Fig.  2D-F).  At  both  constriction  sizes, 
a  dark  circular  arc  corresponding  to  depletion  of  the  fluores- 
cently  labeled  proteins  (henceforth  called  the  depletion  zone) 
on  the  concave  edge  of  the  arc  toward  the  constriction  is  ap¬ 
parent  within  a  fraction  of  a  second  after  the  field  is  turned 
on.  On  the  convex  edge  of  the  arc  away  from  the  constriction, 
we  notice  a  strong  degree  of  preconcentration  of  the  fluores- 
cently  labeled  proteins  (henceforth  called  the  preconcentra¬ 
tion  zone),  especially  close  to  the  sidewall  direction.  It  is  note¬ 
worthy  that  at  both  constriction  sizes,  the  depletion  zone  is 
established  almost  instantaneously  and  does  not  significantly 
evolve  over  time,  while  the  preconcentration  zone  evolves  over 
time.  The  constriction  gap  size  on  the  device  strongly  affects 
the  shape  of  the  arc  and  its  extent  from  the  constriction  (r). 
For  larger  constrictions  (d~140  nm),  the  arc  forms  as  a  semi¬ 
circle,  whereas  for  smaller  constrictions  (~50  nm),  the  arc  is 
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Figure  2.  Time  evolution  of  protein  preconcentration  at  nanoconstriction  devices  of  varying  size  (d  =  50  versus  140  nm)  and  DC  offset  of 
0.3  versus  1.5  V/cm  to  the  NDEP  field  conditions  (200  Vpp/cm,  1  MHz).  (A)-(C)  fluorescence  images  of  50  nm  constriction  devices  under 
1 .5  V/cm  DC  offset  to  NDEP  show  the  instantaneous  formation  of  the  depletion  zone  (~20  p.m)  away  from  the  device  centerline  and  rapid 
emergence  of  the  concentration  zone  close  to  the  sidewall  vector.  (D)-(F)  images  of  140  nm  constriction  devices  under  1 .5  V/cm  DC  offset 
to  NDEP  show  the  instantaneous  formation  of  a  smaller  depletion  zone  (~4  p,m)  away  from  the  device  centerline  and  gradual  emergence 
of  the  concentration  zone  close  to  the  sidewall  vector.(G)  Shows  that  the  preconcentrated  protein  is  immediately  dispersed  upon  turning 
OFF  the  fields.  (H)  Shows  that  a  subcritical  DC  offset  of  0.3  V/cm  is  unable  to  form  a  large  depletion  zone  (r  <1  p.m).  (I)  Shows  that  the 
lower  NDEP  force  at  AC  frequencies  of  approximately  10  MHz,  prevent  formation  of  the  depletion  zone. 
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Figure  3.  Line  scan  fluorescence  intensity  plots  under  NDEP  conditions  with  additional  DC  bias  (~0.3  V/cm  or  ~1 .5  V/cm)  and  varying 
constriction  size  (~50  nm  or  ~140  nm)  along:  (A)  the  centerline  vector  (Y);  (B)  the  sidewall  vector  (Y'). 


extended  along  the  length  (Y-axis)  versus  the  width  (X-axis)  of 
the  channel.  The  extent  of  the  arc  (r),  measured  along  the  de¬ 
vice  centerline  (Y)  was  exponentially  enhanced  from  approxi¬ 
mately  3-4  p.m  for  approximately  140  nm  constriction  gaps  to 
approximately  8-10  p,m  for  100  nm  constriction  gaps  and  ap¬ 
proximately  20-22  p,m  for  approximately  50  nm  constriction 
gap  devices,  as  shown  in  Fig.  S4  of  Supporting  Information. 
At  early  time  points,  protein  preconcentration  is  most  signifi¬ 
cant  along  the  sidewall  direction  (Y'  in  Fig.  1C)  and  over  time 


it  extends  cover  rest  of  the  convex  edge  of  the  arc  facing  the 
inlet,  with  an  almost  unchanging  depletion  region  on  the  con¬ 
cave  edge  of  the  arc  facing  the  constriction.  From  Fig.  2G,  it  is 
apparent  that  upon  removal  of  the  field,  the  preconcentrated 
protein  molecules  are  completely  dispersed  across  the  device, 
demonstrating  a  completely  reversible  preconcentration  pro¬ 
cess,  without  any  discernible  aggregation.  At  a  subcritical  DC 
field  offset  of  0.3  V/cm  to  the  AC  field  under  NDEP  conditions 
(■~200  Vpp/cm  at  1  MFlz),  as  shown  in  Fig.  2H,  the  depletion 
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Figure  4.  Time  evolution  of  the  fluores¬ 
cence  intensity  plots  show  that  an  expo¬ 
nential  level  of  protein  preconcentration 
can  be  achieved  only  at  devices  with  ap¬ 
proximately  50  nm  constriction  gaps  un¬ 
der  NDEP  conditions  with  a  critical  DC 
offset  of  1.5  V/cm.  Legends  from  top  to 
bottom  correspond  to  curves  from  left  to 
right.  The  respective  concentration  lev¬ 
els  for  the  respective  fluorescence  in¬ 
tensities  as  per  the  calibration  curve  are 
indicated. 


zone  has  a  very  small  extent  (<1  |xm)  and  the  preconcentra¬ 
tion  zone  is  wide.  Furthermore,  significant  depletion  and 
preconcentration  zones  are  not  observed  under  AC  fre¬ 
quency  conditions  that  result  in  subcritical  NDEP,  such 
as  at  10  MHz  in  Fig.  21  or  in  the  approximately  30- 
300  kHz  range,  where  only  PDEP  is  apparent  (Fig.  S2, 
Supporting  Information).  Line  scans  of  the  fluorescence  in¬ 
tensity  plots  along  the  centerline  and  sidewall  vectors,  in 
Fig.  3A  and  B,  respectively,  show  that  for  devices  with  ap¬ 
proximately  50  nm  constriction  gaps  under  NDEP  condi¬ 
tions  with  approximately  1.5  V/cm  DC  field,  the  depletion 
zone  is  far  more  extended  along  the  centerline  versus  the 
sidewall  vector,  with  a  far  higher  degree  of  protein  precon¬ 
centration  along  the  sidewall  vector.  On  the  other  hand,  for 
the  same  device  under  NDEP  conditions  plus  subcritical  DC 
(~0.3  V/cm),  the  preconcentration  is  uniform  across  the 
force  balance  arc,  regardless  of  the  centerline  or  sidewall 
direction.  At  larger  constrictions  (~140  nm),  preconcentra¬ 
tion  is  uniform  across  the  arc,  under  NDEP  conditions  plus 
1.5  V/cm  DC  field.  In  summary,  the  results  show  an  expo¬ 


nentially  enhanced  extent  of  protein  depletion  (r)  upon  addi¬ 
tion  of  a  critical  DC  field  to  the  NDEP  conditions  (Fig.  S4, 
Supporting  Information),  especially  at  smaller  constriction 
gaps. 

Figure  4  shows  the  degree  of  protein  preconcentration 
within  devices  of  varying  constriction  sizes  (~50,  100,  or 
140  nm),  AC  frequency  (100  kHz  for  PDEP  and  1  MHz  for 
NDEP),  and  additional  DC  field  offsets  of  approximately  0.3 
or  approximately  1.5  V/cm.  The  results  are  plotted  as  maxi¬ 
mum  fluorescence  intensity  (localized  on  the  convex  edge  of 
the  arc  in  proximity  to  the  constriction  tip)  versus  trapping 
time.  The  approximate  local  concentration  level  at  critical  flu¬ 
orescence  intensity  levels  as  per  the  concentration  calibration 
curve  is  indicated  for  reference.  From  this  figure,  it  is  clear 
that  under  NDEP  conditions  with  subcritical  DC  field  (~0.3 
V/cm),  the  degree  of  protein  preconcentration  is  not  substan¬ 
tial,  even  for  devices  with  approximately  50  nm  constriction 
gaps.  However,  the  addition  of  approximately  1.5  V/cm  DC 
field  to  the  NDEP  conditions  causes  a  rapid  rise  in  the  de¬ 
gree  of  protein  preconcentration,  far  above  what  is  obtained 
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Figure  5.  Simulations  of  the  net  force 
(F net  =  Fdep  -  Fep  +  FEo)  along  the  center- 
line  (V)  and  sidewall  ( Y ')  vectors  (as  per 
Figure  1C)  show  a  sharp  drop-off  in  the 
force  along  the  sidewall  versus  the  cen¬ 
terline  vector.  The  inset  shows  the  pro¬ 
file  across  the  device  of  the  electrophore¬ 
sis  (blue  arrow)  versus  NDEP  (red)  and 
electro-osmosis  (green)  forces. 
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under  PDEP  conditions  inside  the  constriction  region.  The 
rise  in  degree  of  preconcentration  is  steep  for  50  nm  con¬ 
striction  devices,  where  the  detector  is  saturated  within  2  s, 
based  on  multiple  data  points  acquired  at  30  frames  per  sec¬ 
ond.  In  summary,  these  results  indicate  the  enhanced  extent 
of  the  protein  depletion  zone  under  conditions  of  a  char¬ 
acteristic  AC  field  (200  Vpp/cm  at  1  MHz)  and  a  critical  DC 
offset  of  approximately  1.5  V /cm  (as  in  Fig.  2)  on  devices  with 
approximately  50-nm  constriction  gaps,  can  be  directly  corre¬ 
lated  to  rapid  protein  preconcentration  (Fig.  4). 

4  Discussion 
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Figure  6.  Simulation  of  potential  energy  profiles  to  explain  for¬ 
mation  of  the  depletion  and  preconcentration  zones.  (A)  Under 
NDEP  conditions  (200  Vpp/cm,  1  MHz)  at  50  nm  constriction  gap 
devices,  a  barrier  of  substantially  high  potential  energy  can  be 
formed.  (B)  Upon  the  addition  of  a  DC  offset  of  1.5  V/cm  to  the 
NDEP  conditions,  the  barrier  is  tilted  toward  the  inlet  to  form  po¬ 
tential  energy  wells  for  protein  preconcentration.  (C)  The  sharp 
drop-off  in  potential  energy  away  from  the  constriction  along  the 
sidewall  versus  the  centerline  vector  is  especially  apparent  within 
devices  of  smaller  (~50  nm)  versus  larger  (~500  nm)  constriction 
gaps.  (D)  Under  a  critical  DC  offset  under  NDEP  conditions,  the 
potential  profile  is  tilted  to  cause  deeper  and  narrower  wells  for 
protein  preconcentration. 


The  instantaneous  formation  of  the  protein  depletion  zone  of 
enhanced  extent  (~20  p,m)  in  the  vicinity  of  the  constriction 
of  the  approximately  50-nm  gap  devices  and  its  direct  correla¬ 
tion  to  the  steeply  rising  degree  of  protein  preconcentration 
is  qualitatively  explained  in  this  section  through  applying 
device-level  simulations  to  describe  the  electrokinetic  force 
balance.  In  comparison,  the  depletion  region  reported  within 
prior  work  [4-6]  due  to  concentration  polarization  or  electri¬ 
cal  double  layer  overlap  at  the  micro-to-nanofluidic  interface 
evolves  over  time  and  hence,  results  in  a  slower  degree  of 
preconcentration.  Furthermore,  we  seek  to  explain  the  dis¬ 
parities  in  protein  preconcentration  along  the  sidewall  versus 
centerline  directions  that  are  observed  only  within  the  smaller 
constriction  devices  (~50  nm  gap)  under  NDEP  conditions 
(200  Vpp/cm,  1  MHz)  with  a  critical  DC  field  (1.5  VDC/cm),  as 
in  Fig.  3.  Finally,  it  is  of  interest  to  correlate  these  preconcen¬ 
tration  disparities  to  the  unexpectedly  enhanced  extent  of  the 
depletion  zone  (~20  p,m  in  Fig.  2A-C)  within  these  devices, 
given  that  the  influence  of  the  NDEP  force  should  drop  off 
within  an  extent  (r)  of  approximately  1  p,m  away  from  the  con¬ 
striction  edge  and  the  total  extent  of  the  constriction  gradient 
along  the  centerline  vector  is  no  more  than  5  |xm. 

The  inset  of  Fig.  5  shows  a  simulation  of  the  direction 
and  magnitude  (length  of  the  arrows)  of  the  electrophoresis 
(blue),  DEP  (red),  and  electro-osmosis  (green)  forces  across 
the  device.  It  is  readily  apparent  that  similar  to  the  schematic 
analysis  presented  in  Fig.  1C,  the  electrophoresis  and  electro¬ 
osmosis  forces  oppose  each  other  at  all  points  within  the 
device,  and  both  of  these  forces  rise  continuously  with  the 
electric  field  due  to  gradient  neighboring  the  constriction,  to 
reach  a  maximum  value  at  the  constriction  tip.  FEP  exceeds 
the  Feo  at  all  points  within  the  device,  with  Fndep  making 
up  the  difference  in  the  region  proximal  to  the  constriction 
to  cause  preconcentration  on  the  side  facing  the  inlet  as  per 
the  force  balance  condition:  FEP  =  FNdep  +  Feo.  Hence,  in 
the  absence  of  a  critical  level  of  Fndep,  no  significant  protein 
preconcentration  would  occur,  as  apparent  from  the  exper¬ 
imental  results  in  Fig.  21,  where  the  preconcentration  zone 
is  insignificant  under  NDEP  conditions  at  10  MHz,  due  to 
lower  Fndep  in  comparison  to  that  at  1  MHz.  Furthermore, 
while  FEP  and  FEO  scale  with  the  field,  FNDEP  scales  as  VE2. 
Hence,  based  on  field  simulations,  FNDEP  rises  more  sharply 
in  the  vicinity  of  the  50-nm  constriction  gap  to  a  far  higher 
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Figure  7.  Schematics  for  formation  of  deple¬ 
tion  and  preconcentration  zones:  (A)  NDEP 
conditions  with  a  critical  DC  field  offset  at 
50  nm  constriction  gaps  causes  large  deple¬ 
tion  zones  and  more  localized  preconcentra¬ 
tion  zones  due  to  the  large  potential  energy 
disparities  along  centerline  versus  sidewall  di¬ 
rection.  The  respective  potential  energy  dis¬ 
parities  are  relatively  lower  within  devices  of: 
(B)  larger  constriction  gaps,  due  to  broader  po¬ 
tential  wells;  and  (C)  under  conditions  of  sub- 
critical  DC,  due  to  shallower  potential  wells, 
thereby  causing  depletions  zones  of  lower  ex¬ 
tent  and  less  localized  preconcentration  zones. 


value  (maximum  of  ~10~11  N)  than  either  FEP  (maximum  of 
~10~12  N)  or  FE0  (maximum  of  ~10~13  N),  but  FNDEP  also 
drops  off  more  sharply  away  from  the  constriction  tip.  To 
better  understand  the  effect  of  this  on  the  electrokinetic  force 
balance  across  the  device,  the  net  electrokinetic  force  (Fnet, 
as  in  Eq.  (6)  is  plotted  in  Fig  5.,  along  the  centerline  versus 
sidewall  directions  for  large  (~500  nm)  and  small  (~50  nm) 
constriction  gaps.  Approaching  from  the  inlet  (right  side  of 
the  plot)  toward  the  constriction  (left  side  of  the  plot),  the  Fnet 
rises  more  sharply  along  the  sidewall  versus  the  centerline 
direction  due  to  the  gradient  resulting  from  the  constriction, 
and  this  rise  is  sharper  for  the  smaller  constriction  gap  sizes. 
Based  on  the  length  of  the  arrows  within  the  inset  of  Fig.  5, 
representing  magnitude  of  the  respective  forces,  it  is  clear  that 
at  an  extent  of  r  approximately  1  p,m  from  the  constriction, 
while  Fndep  almost  equals  FEP  along  the  centerline  direction; 
it  is  insignificant  in  comparison  FEP  along  the  sidewall  direc¬ 
tion.  Hence,  approximately  1  p,m  away  from  the  constriction, 
the  steeper  drop  in  Fnet  along  the  sidewall  versus  centerline 
direction  is  directly  attributable  to  the  sharper  drop  off  in 
Fndep  along  the  sidewall  versus  centerline  direction. 

Next,  we  relate  these  force  profiles  to  the  scalar  potential 
energy  across  the  device  to  understand  and  optimize  device 
designs  for  enhanced  preconcentration,  similar  to  prior  work 
utilizing  this  methodology  [30].  Specifically,  we  focus  on  ap¬ 
plying  this  approach  to  qualitatively  understand  the  shape 
and  extent  of  the  depletion  region,  as  well  as  the  location  of 
preconcentration  regions  within  the  device.  In  Fig.  6A  and 
B,  the  potential  energy  profile  due  to  only  FNDEP  (UDEp)  and 
Fnet  ( U)  is  plotted  in  KBT  units,  across  from  the  entire  device, 
from  the  inlet  to  the  constriction  to  the  outlet.  Under  AC 
field  conditions  (only  NDEP),  a  potential  energy  barrier  due 
to  Fndep  develops  in  close  proximity  of  the  constriction,  which 
prevents  protein  particles  from  passing  through  the  constric¬ 
tion.  Upon  addition  of  a  critical  level  of  DC  field  offset,  the 
potential  energy  profile  is  tilted  to  result  in  “wells”  for  more 
localized  protein  trapping.  The  influence  of  the  resulting 
wells  across  the  centerline  and  sidewall  directions  is  explored 
for  devices  of  varying  constriction  gap  size  and  DC  bias  in 
Fig.  6C  and  D.  As  the  constriction  is  approached  onward 
from  the  inlet,  disparities  emerge  in  the  potential  energy  bar¬ 
rier  along  the  sidewall  versus  centerline  directions  in  Fig.  6C. 
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These  disparities  are  enhanced  at  smaller  constriction  gap 
size  due  to  the  substantially  higher  Fndep  within  smaller  gaps, 
which  manifests  as  a  higher  Fnet  or  negative  slope  of  the  po¬ 
tential  energy  profile  along  the  sidewall  direction.  Figure  6D 
shows  that  a  greater  DC  field  offset  to  tilt  the  potential  energy 
profile  can  cause  deeper  potential  wells  for  protein  precon¬ 
centration  from  the  inlet,  as  well  as  from  regions  of  higher 
energy  across  the  entire  device,  including  those  points  along 
the  centerline  direction  that  have  shallower  and  wider  poten¬ 
tial  wells  (due  to  smaller  changes  in  Fnet  along  the  centerline 
versus  sidewall  direction,  as  per  Fig.  5).  Hence,  within  de¬ 
vices  of  smaller  constriction  gaps,  at  a  critical  DC  field  offset, 
the  sidewall  direction  of  the  device  has  deepest  and  narrowest 
potential  wells.  This  results  in  depletion  of  protein  molecules 
from  all  directions  other  than  the  sidewall.  The  depletion  is 
greatest  along  the  centerline  direction,  which  has  the  largest 
potential  energy  disparities  from  that  of  the  sidewall  direction. 
As  per  Fig.  7A,  at  a  critical  DC  field  offset  under  NDEP  condi¬ 
tions,  the  ensuing  depletion  can  cause  net  preconcentration 
along  the  sidewall  direction  within  a  few  microns  from  the 
constriction  edge.  For  devices  with  larger  constriction  gaps, 
the  potential  energy  disparities  between  the  sidewall  and  cen¬ 
terline  directions  are  lower,  as  in  Fig.  6C,  thereby  resulting 
in  lower  preconcentration  at  the  wider  potential  wells,  as  in 
Fig.  6D.  This  explains  the  experimental  data  indicating  more 
uniform  preconcentration  along  the  centerline  (in  Fig.  3A) 
versus  the  sidewall  direction  (in  Fig.  3B)  for  the  140-nm  con¬ 
strictions,  as  opposed  to  the  larger  differences  along  these  di¬ 
rections  for  50-nm  constriction  devices.  It  also  causes  a  lower 
degree  of  depletion  at  larger  constriction  gaps  (Fig.  7  A  versus 
Fig.  7B)  and  explains  the  lower  extent  of  depletion  zone  in 
the  experimental  data  of  Fig.  2D  (140-nm  constriction)  ver¬ 
sus  Fig.  2A  (50-nm  constriction).  At  subcritical  DC  fields,  the 
potential  wells  are  not  deep  enough  to  cause  substantial  pre¬ 
concentration  in  comparison  to  the  case  where  critical  DC 
conditions  are  applied.  This  causes  a  lower  degree  of  deple¬ 
tion  and  preconcentration  as  per  the  schematics  in  Fig.  7C. 
This  is  consistent  with  the  experimental  data  in  Fig.  2H,  as 
well  as  more  uniform  preconcentration  along  the  centerline 
versus  sidewall  (Fig.  3A  and  B)  at  conditions  of  subcritical 
DC  (0.3  V/cm)  versus  critical  DC  field  (1.5  V/cm).  Hence, 
sub  100-nm  constriction  gaps  enhance  the  energy  barrier  at 
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the  constriction  and  create  more  disparities  at  the  centerline 
versus  sidewall  direction  of  the  device  to  enable  narrow  po¬ 
tential  wells  along  the  sidewall  direction  for  rapid  and  steep 
degrees  of  protein  preconcentration.  A  critical  DC  field  offset 
to  sufficiently  tilt  the  potential  profile  coupled  to  significantly 
small  constriction  gaps  is  required  to  cause  deep  and  narrow 
potential  wells  for  protein  preconcentration.  Finally,  a  crit¬ 
ical  level  of  Fndep,  which  only  occurs  at  the  characteristic 
frequency  of  1  MHz  for  streptavidin  protein  molecules,  is  re¬ 
quired  to  localize  the  preconcentration  region  in  the  vicinity 
of  the  nano-  constriction. 

In  conclusion,  a  methodology  utilizing  a  balance  of  elec- 
trokinetic  forces  at  nano-constriction  devices  is  described  for 
enabling  steeply  rising  levels  of  preconcentration  of  strepta¬ 
vidin  protein  molecules  within  physiological  media  based  on 
the  following  three  requirements.  First,  a  characteristic  AC 
frequency  of  1  MHz  is  required  to  enable  trapping  of  the 
protein  molecules  under  negative  dielectrophoresis  (NDEP) 
conditions  using  a  200  Vpp/cm  field.  Second,  devices  with 
approximately  50-nm  constriction  gaps  are  required  to  cause 
narrow  potential  energy  wells  along  the  sidewall  direction 
of  the  device.  Third,  a  critical  level  of  DC  field  offset  to  the 
AC  field  (experimentally  determined  to  be  1.5  V/cm  for  ap¬ 
proximately  50-nm  constriction  gap  devices)  is  required  to 
sufficiently  tilt  the  potential  profile  to  cause  deep  potential 
energy  wells  for  protein  preconcentration.  Under  these  con¬ 
ditions,  an  elliptical-shaped  depletion  zone  of  larger  extent 
along  the  device  centerline  length  axis  (~20  pan  for  ~50- 
nm  constriction  gap  devices)  than  other  directions  forms  in¬ 
stantaneously  around  the  constrictions  to  result  in  rapidly 
rising  levels  of  protein  preconcentration  along  the  sidewall 
direction  from  the  constriction.  While  protein  preconcen¬ 
tration  can  occur  at  NDEP  conditions  within  devices  con¬ 
taining  larger  constrictions  (up  to  ~140  nm)  and  under 
conditions  of  lower  DC  field  offset,  the  smaller  disparities 
in  the  potential  energy  profile  along  the  sidewall  versus 
centerline  directions  cause  depletion  zones  of  smaller  ex¬ 
tent  to  result  in  a  smaller  rise  in  the  degree  of  preconcen¬ 
tration  over  time.  These  protein  preconcentration  method¬ 
ologies  may  be  applied  toward  biomarker  discovery,  pro¬ 
tein  crystallization,  and  rare  target  sensing  for  early  disease 
diagnostics. 
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Scaling  analysis  of  constriction-based  dielectrophoresis  devices  for  trapping  nanoscale  bio¬ 
particles  in  physiological  media  of  high-conductivity 
ABSTRACT 

Selective  trapping  and  pre-concentration  of  nanoscale  bio-particles  (size  <100  rnn)  is  significant 
for  the  separation  and  high-sensitivity  detection  of  biomarkers.  Dielectrophoresis  is  capable  of 
highly  selective  trapping  of  bio-particles  based  on  their  characteristic  frequency  response. 
However,  the  trapping  forces  fall  steeply  with  particle  size,  especially  within  physiological 
media  of  high-conductivity  where  the  trapping  can  be  dissipated  by  electrothermal  flow  due  to 
localized  Joule  heating.  Herein,  we  investigate  scaling  of  the  device  within  an  electrodeless 
dielectrophoresis  geometry  through  the  application  of  highly  constricted  microchannels  coupled 
to  the  reduction  of  channel  depth  to  explore  the  net  balance  of  dielectrophoretic  trapping  force 
versus  electro  thermal  drag  force  on  bio-particles.  While  higher  degrees  of  constriction  enable 
dielectrophoretic  trapping  of  successively  smaller  bio-particles  within  a  short  time,  the 
electrothennal  flow  due  to  enhanced  Joule  heating  within  media  of  high  conductivity  can  cause  a 
significant  dissipation  of  bio-particle  trapping.  The  dissipative  effect  of  the  electrothennal  drag 
force  can  be  reduced  through  lowering  the  depth  of  the  highly  constricted  channels  to  sub-micron 
sizes,  which  substantially  reduces  the  degree  of  Joule  heating,  thereby  enhancing  the  range  of 
voltages  and  media  conductivities  that  can  be  applied  towards  rapid  dielectrophoretic  pre¬ 
concentration  of  silica  nanoparticles  (~50  rnn)  and  streptavidin  protein  biomolecules  (~5  nm). 

We  envision  the  application  of  these  methodologies  towards  nanofabrication,  optofluidics, 
biomarker  discovery,  and  early  disease  diagnostics. 

Keywords:  Dielectrophoresis,  Joule  heating,  Electrothennal  flow,  Nanoanalysis,  Proteins 


Scaling  analysis  of  constriction-based  dielectrophoresis  devices  for  trapping  nanoscale  bio¬ 
particles  in  physiological  media  of  high-conductivity 
I.  Introduction 

High-sensitivity  detection  of  rare  numbers  of  biomarkers  requires  methodologies  for 
selective  pre-concentration  of  the  biomarker  within  physiologically  relevant  media  in  the 
proximity  of  the  sensor  [1],  In  comparison  to  chemical  affinity  methods  based  on  antibody 
depletion  [2],  electrokinetic  methods  within  micro-  or  nanofluidic  devices,  can  achieve  far  higher 
degrees  of  pre-concentration  due  to  the  large  volume  reduction  [3].  Dielectrophoresis  (DEP) 
enables  highly  selective  trapping  of  bio-particles  based  on  the  characteristic  frequency  response 
of  the  dielectric  permittivity  of  the  bio-particle  versus  that  of  the  medium  [4],  [5],  and  it  has  been 
extensively  applied  towards  sorting  of  somewhat  similar  sized  biological  cells  with  differing 
dielectric  frequency  response  [6].  However,  its  application  towards  trapping  5-50  nm  sized 
biomarkers  requires  structures  such  as  insulator  posts  or  dielectric  constrictions  to  enhance  the 
local  field  using  an  electrodeless  dielectrophoresis  scheme,  to  offset  the  steep  fall  in  DEP 
trapping  forces  with  particle  size  [7],  [8],  [9],  [10],  as  demonstrated  for  trapping  proteins  [11], 
[12],  [13],  [14],  [15]  and  fragments  of  ss-DNA  [16],  [17].  Biomarkers  need  to  be  trapped  under 
room  temperature  conditions,  within  physiological  media  of  neutral  pH  and  high  conductivity,  to 
maintain  their  conformation  and  bio-functionality.  Hence,  while  dielectric  constrictions  enhance 
dielectrophoretic  trapping  forces,  they  also  cause  a  significant  degree  of  Joule  heating  within 
media  of  high  conductivity  (>1  S/m)  due  to  enhanced  current  densities  at  the  high- field  points 
[18],  [19].  The  resulting  localized  temperature  rise  can  adversely  affect  biomarker  functionality 
(such  as  causing  protein  denaturation)  and  the  ensuing  temperature  gradients  can  cause 
electrothennal  flow  that  dissipates  dielectrophoretic  trapping. 


The  generation  term  for  electrothermal  flow:  amE2,  where  am  is  media  conductivity  and  E  is 
the  local  electric  field,  acts  over  a  significantly  long-range,  due  to  its  indirect  action  on  the 
particles  through  a  drag  force.  The  trapping  force  due  to  DEP,  on  the  other  hand,  is  highly 
localized  at  the  constriction  due  to  its  dependence  on  VE2.  Hence,  within  the  constriction  device 
geometry,  electrothermal  flow  acts  to  dissipate  trapping  under  positive  DEP,  due  to  their 
opposing  directions  [9];  while  under  negative  DEP  conditions,  it  causes  spreading  of  the  trapping 
zone  over  a  diffuse  region,  thereby  reducing  the  degree  of  pre-concentration  [11].  While  the  DEP 
trapping  force  scales  as:  VE2,  the  low-frequency  (<  1  MHz)  Coulombic  component  of 
electrothennal  flow  scales  as:  (VT.E)  E;  where  VT  is  the  temperature  gradient  [20],  Hence,  a 
scaling  analysis  of  the  DEP  force  versus  electrothermal  drag  force  for  varying  trapping 
conditions  (applied  field  and  media  conductivity)  and  device  conditions  (constriction  ratio  and 
channel  depth)  can  lead  to  nanofluidic  device  designs  that  are  optimized  to  enhance  trapping 
over  dissipation.  In  this  current  work,  we  study  how  scaling  of  the  device  geometry,  media 
conductivity  and  applied  voltage  affect  DEP  forces,  localized  temperatures  and  the  balance  of 
DEP  trapping  (FDep)  versus  electrothennal  drag  force  (FEt)  on  nanoscale  bio-particles. 

II.  METHODOLOGY 

A.  Device  Geometry:  The  device  geometry  of  a  dielectric  constriction  is  shown  in  Figure  1.  For 
experiments,  the  constriction  channels  were  either  patterned  on  PDMS  or  etched  into  a  quartz 
substrate  and  then  bonded  onto  cover  slip  glass  for  imaging  the  DEP  translation  of  fluorescently 
labeled  streptavidin  protein  (52.8  kDa,  Molecular  Probes,  Engene,  OR)  or  -50  nm  silica 
nanoparticles  (Corpuscular  Inc.,  Cold  Spring,  NY)  using  an  inverted  microscope  (Zeiss  Observer 
Al).  The  depth  of  the  channel  is  indicated  for  each  of  the  results. 


B.  Device  Modeling:  We  focus  herein  on  the  balance  of  dielectrophoresis  trapping  force  versus 
electrothennal  flow  for  microfluidic  channel  constriction  ratios  varied  from  50x  (500  pm 
reduced  to  10  pm)  to  5000x  (500  pm  reduced  to  0.1  pm);  media  conductivity  varied  from:  0.1-10 
S/m;  and  the  applied  field  varied  from:  100-400  V/cm.  Finite  element  methods,  using  the  ESI- 
CFD  ACE+  software  (ESI  CFD,  North  America,  Huntsville,  Al)  were  used  for  the  computations. 

Dielectrophoresis:  The  constriction  micro-channel  device  works  on  the  principle  that  the  highly 
insulating  constriction  material  enhances  the  applied  field  at  the  constriction  tip.  The  average 
steady-state  dielectrophoretic  force,  Fdep,  on  a  homogeneous  spherical  particle  with  electrical 
pennittivity  sP;  conductivity  ap  and  radius  a  suspended  in  a  fluid  with  electrical  pennittivity  sm 
and  conductivity  cm  in  a  non-uniform  AC  electric  field  E  is  given  by  [4]: 

2 

..(Eq.  1) 

Here,  KCm  is  the  Clausisus-Mossotti  factor  which  depends  on  the  frequency  (to)  of  the  applied 
electric  field  as:  s*  =  s  +  a/jco.  The  DEP  force  is  experienced  by  the  polarized  particles  only 
under  non-uniform  electric  fields;  hence,  it  is  proportional  to  the  gradient  (or  spatial  non¬ 
uniformity)  of  the  square  of  the  electric  field.  For  computational  simplicity,  we  herein  assume 
that  Kcm  is  ~1,  even  though  the  value  can  range  from  0  to  1,  for  experimental  positive  DEP 
conditions. 

Joule  heating  and  temperature  profile:  Electrical  conductivity  of  the  fluidic  media  in  the  channel 
under  the  applied  field  during  dielectrophoresis,  leads  to  current  flow  and  hence  Joule  heating, 
thereby  causing  a  temperature  rise.  Since  the  electric  field  in  the  channel  is  not  uniform,  the 
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heating  is  also  non-uniform.  At  steady  state,  heat  generated  in  the  system  is  equal  to  the  heat 
dissipated.  Thus  the  temperature  profde  can  be  computed  as  [5]: 


kV2T  + 


=  0 


...(Eq.  2) 


Where  k  is  the  thennal  conductivity  of  the  fluid,  T  the  temperature  and  o|E|2  is  the  Joule  heating 
tenn.  Here  the  convection  term  of  heat  transfer  is  neglected  owing  to  the  micro-scale  geometry 
of  the  device  where  thermal  conduction  is  main  mechanism  of  heat  transfer. 


Electrothermal  force:  Electrothennal  flow  is  the  body  force  (Fet)  experienced  by  the  fluid  in 
presence  of  an  applied  electric  field  due  to  localized  anisotropic  heating  in  the  channel.  Gradients 
in  temperature  in  the  channel  give  rise  to  gradients  in  the  electrical  properties  of  the  media, 
namely  conductivity  and  dielectric  permittivity.  For  an  incompressible  fluid,  this  is  [20]: 


Fet  ~  PF  ^ 


Vf 


...(Eq.  3) 


Here,  p  is  the  charge  density  (as  a  result  of  temperature  gradients),  s  is  the  dielectric  permittivity 
and  E  is  the  applied  electric  field.  The  first  term  on  the  right  hand  side  of  the  equation  is  the 
Coulomb  force  (Fc)  while  the  second  term  is  the  dielectric  force  (Fd).  Since,  the  origin  of  these 
forces  is  in  the  temperature  gradients,  they  can  be  simplified  and  expressed  as  the  time  averaged 
force  per  unit  volume  in  terms  of  temperature  gradients  as  [5]: 
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...(Eq.  4) 


As  apparent  from  the  expression  above,  the  electrothermal  forces  are  frequency  dependent.  At 
low  frequency,  the  Coulomb  forces  (Fc)  are  dominant  while  at  high  frequencies  the  Dielectric 


forces  (Fd)  become  dominant  [20].  We  focus  on  low  frequency  AC  fields  (<  1  MHz)  within  this 
work;  hence,  the  dominant  electrothermal  force  is  the  Coulomb  force  (Fc).  As  per  [21],  the 

temperature  dependence  of  the  conductivity  and  pennittivity  are  given  by  2%  and 

~zr  —  —0.4%  .  Hence,  since:  -zr  ~  ~^r  >  0,  the  force  direction  can  be  estimated  based  on  the 

£  dT  ’  a  dT  as  dT  ’ 

sign  of VT  •  E.  When  the  temperature  gradient  is  parallel  to  the  electric  field,  the  Coulomb  force 
is  directed  against  the  electric  field  and  vice  versa. 

Electrothermal  Flow  Profile:  Since  the  body  force  due  to  electrothermal  flow  acts  to  move  the 
bulk  fluid,  it  can  be  included  within  the  Navier-Stokes  equation.  Thus,  given  the  body  force  and 
assuming  an  incompressible  fluid  with  low  Reynolds  number,  the  steady  state  Navier-Stokes 
equation  can  be  written  as: 

r]V2v  —  Vp  +  (F)  =  0  ...(Eq.  5) 

— > 

Here,  77  is  the  viscosity,  p  is  the  pressure  and  ( F )  is  the  average  volumetric  force  on  the  fluid.  In 

this  case  (F)  can  be  approximated  as  [22]:  (F)  =  (Fc  +  FD).  As  the  pressure  in  the  fluid  is 

constant  and  at  low  frequencies,  electrothennal  force  is  dominated  by  the  Coulomb  force,  the 
equation  is  reduced  to 

77V2v  +  <Fc>  =  0  ...(Eq.  6) 

Balance  of  forces  on  the  bio-particles:  Finally,  we  examined  the  balance  of  FET  versus  FDEP 
forces  on  the  net  mobility  for  the  bio-particles  for  varying  device  conditions  (constriction  ratio 
and  channel  depth)  and  field  conditions  (applied  field  and  media  conductivity).  The  bulk  fluid 
flow  acts  on  the  bio-particle  in  the  fonn  of  drag  force,  thereby  affecting  the  particle  trajectory. 
Particle  tenninal  velocity  in  the  fluid  can  be  given  by: 


...(Eq.  7) 
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Here,  vpis  the  net  particle  velocity,  Vf  is  the  fluid  velocity  due  to  electrothennal  flow,  Fp  is  the 

— >  — > 

force  acting  on  the  particle,  which  in  this  case  is  the  dielectrophoretic  force  ( Fp  —  FDEP  )  and  y  is 
the  drag  coefficient  for  the  particle.  Assuming  a  spherical  particle  of  radius, a,  the  drag 
coefficient  is  given  by :y  =  6717712...  (Eq.  8) 

From  this  analysis,  net  force  on  the  particle  can  be  computed  as: 

FpTotal  —  Yvf  +  Fdep 

. ..(Eq.  4 

III.  Results 

TIT. I.  Enhancing  DEP  trapping  through  scaling  the  constriction  ratio:  Based  on  the  cubic 
dependence  of  the  DEP  force  on  the  hydrodynamic  radius  of  the  particle  (a),  as  per  (Eq.  1),  the 
DEP  force  decreases  cubically  with  particle  size.  As  per  the  device  structure  under  the  DEP  field 
in  Figure  1,  an  insulator  constriction  enhances  the  conduction  current  (at  low  frequencies)  and 
displacement  current  (at  high  frequencies)  within  the  constricted  gap  region  of  the  device, 
thereby  enhancing  the  density  of  field  lines  in  the  vicinity  of  the  constriction  edge  to  enable  the 
localized  enhancement  of  DEP  trapping  forces  due  to  higher  field  (dependent  on  constriction 
ratio,  R)  and  its  gradient  (dependent  on  sharpness  of  the  constriction).  We  consider  a  field  of  350 
V/cm  within  channels  containing  constrictions  with  ratios  (R)  of  500x  (500pm  to  1pm)  and 
5000x  (500pm  to  0.1pm)  for  trapping  bio-particles  of  sizes  50nm  and  5nm,  within  media  of  1 
S/m  conductivity.  Comparing  trapping  forces  to  the  Brownian  motion  on  particles,  the  minimum 
DEP  force  to  trap  50  nrn  and  5  nm  particles  can  be  determined  as  ~1  fN  and  ~3  fN,  respectively, 


as  per  Appendix  A,  whereas  a  significant  degree  of  pre-concentration  within  minutes  would 
require  forces  in  the  pN  range  [Soory].  From  °  (a)  and  (c),  it  is  clear 

that  for  50  nm  bio-particles,  ~pN  range  DEP  force  distributions  can  be  attained  over  a 
considerable  extent  of  the  constriction  region  for  devices  with  500x  and  5000x  fold  constrictions. 
However,  upon  a  one  order  of  magnitude  reduction  in  bio-particle  radius  (from  a  =  50  nm  to  5 
nm),  the  magnitude  of  the  DEP  force  is  lowered  by  3  orders  of  magnitude.  As  can  be  seen  from 
Ih^!  £% 'f'  5!j 0 M  }jgi  o  (b),  the  regions  of  device  with  DEP  force  greater  than  ~pN  are 
almost  negligible  for  devices  with  500x  constrictions,  whereas  with  5000x  constrictions  in 
Figure  2(d),  this  region  with  >  pN  DEP  trapping  forces  is  restricted  to  the  constriction  edges. 
Since  bio-particles  under  the  effect  of  DEP  force  move  with  a  velocity  given  by:  v  —  F /6nr)a, 
this  reduction  in  force  magnitude  will  in  turn  cause  a  reduction  in  particle  trapping  velocity  and  a 
proportionate  increase  in  particle  trapping  time.  Based  on  the  magnitudes  of  the  velocity  profiles 
of  Figure  2(e)-(h),  we  can  estimate  that  instead  of  trapping  times  on  the  order  of  a  minute  for  50 
nm  bio-particles  using  500x  constrictions  [16],  the  same  constriction  ratio  would  likely  take 
100’s  of  minutes  to  achieve  a  similar  degree  of  pre-concentration  with  5nm  bio-particles.  This 
trend  is  substantiated  by  our  experimental  observations  of  streptavidin  protein  pre-concentration 
in  Figure  3,  showing  that  while  a  1000-fold  degree  of  pre-concentration  can  be  obtained  within 
~20  seconds  using  9000x  constrictions  (constriction  from  30  pm  to  50  nm  in  X  and  Y  directions, 
and  from  3  pm  to  200  nm  in  Z  direction,  as  described  previously  [12]),  the  same  degree  of  pre¬ 
concentration  would  likely  take  tens  of  minutes  using  3000x  constrictions  (constriction  from  30 
pm  to  150  nm  in  X  and  Y  directions,  and  from  3  pm  to  200  nm  in  Z  direction),  and  no 
perceptible  pre-concentration  can  be  observed  with  lOOOx  constrictions. 


III.2.  Constraints  on  scaling  imposed  by  Joule  heating:  The  chief  constraint  to  scaling  down 
the  constriction  gap  (higher  R)  arises  from  Joule  heating  in  the  constriction  region.  In  order  to 
maintain  the  functionality  of  the  bio-particle,  such  as  conformations  required  for  its  participation 
within  biochemical  recognition  events,  the  media  conductivity  (om)  needs  to  be  in  the  range  of  1- 
10  S/m  (-50-500  mM  salt  levels).  The  application  of  DEP  fields  within  media  of  high- 
conductivity  causes  an  enhanced  level  of  current  flow,  especially  within  devices  containing 
dielectric  constrictions  to  significantly  enhance  the  electric  field  (E).  This  causes  Joule  heating 
[23],  which  is  proportional  to:  am\E\2.  Furthermore,  the  non-unifonnity  of  the  field  in  the 
constriction  region  versus  elsewhere  in  the  channel  causes  the  Joule  heating  to  be  non-uniform. 
Additionally,  since  the  floor  and  ceiling  of  the  micro-channel  act  as  the  predominant  heat  sinks 
for  the  generated  heat  in  the  fluid  through  Joule  heating,  temperature  gradients  (VT)  are  formed. 
A  hot  spot  develops  midway  through  the  channel  depth  (5  pm),  as  shown  in  Figure  4a,  where 
local  temperatures  can  rise  to  43.4°C  (316.4K)  for  a  5000x  constriction  (500  to  0.1pm 
constriction)  under  350  V/cm  fields  within  media  of  1  S/m  conductivity.  In  Figure  4b  we  show 
evidence  for  the  destruction  of  PDMS  constriction  tips  due  to  significant  levels  of  Joule  heating, 
especially  in  the  central  region  of  channels  of  -  5  pm  depth.  These  significant  levels  of  Joule 
heating  can  cause  the  denaturation  of  proteins,  and  the  resulting  temperature  gradients  drive 
electro  thermal  flow  in  the  bulk  fluid.  As  discussed  previously,  such  electrothermal  flow  can 
dissipate  trapping  under  positive  DEP  conditions  [9]  and  reduce  the  degree  of  pre-concentration 
under  negative  DEP  conditions  [11]. 

To  understand  the  limits  imposed  by  Joule  heating  on  the  range  of  applied  voltages  and 
media  conductivity  that  can  be  applied  towards  dielectrophoretic  trapping  within  devices 
containing  highly  scaled  constriction  ratios  (large  R),  we  present  temperature  simulations  of  the 


hotspot  region  under  these  conditions.  It  is  apparent  here  that  the  Joule  heating  causes  a  non¬ 
linear  rise  in  the  temperature  with  media  conductivity  (Figure  4c)  and  with  applied  field  (Figure 
4d).  With  500x  constriction  devices  of  5  pm  depth,  the  temperature  rises  to  just  a  few  degrees 
above  the  ambient  (300  K)  for  the  voltages  studied  herein,  as  per  Figure  4d.  However,  with 
5000x  constriction  devices  of  5  pm  depth,  the  temperature  rises  sharply  to  several  tens  of 
degrees  above  ambient  at  fields  >  250  V/cm  within  media  of  om  =  1  S/m  (Figure  4d);  and  the 
temperature  can  rise  to  aqueous  boiling  point  levels  at  am  =  10  S/m,  even  at  100  V/cm  fields 
(Figure  4c).  With  decreasing  channel  depth  (t),  the  generation  term  due  to  Joule  heating  can  be 
significantly  reduced  due  to  the  dependence  of  temperature  at  the  hotspot  (Tmax)  on  the  square  of 
channel  depth  (t2),  based  on  a  parallel  plate  capacitor  model  (Appendix  B).  Hence,  a  significant 
lowering  of  the  temperature  at  the  hotspot  is  apparent  in  Figure  4c  and  d. 

III.3.  Balance  of  dielectrophoretic  trapping  force  versus  electrothermal  drag  force: 

Temperature  gradients  within  the  constriction  device  can  cause  a  significant  level  of 
electro  thermal  flow  on  the  fluid,  which  can  influence  the  overall  particle  trajectory  and  lead  to 
increased  drag  force  on  the  particles.  In  the  vicinity  of  the  hotspot  region  within  the  constriction 
(Figure  4a),  the  direction  of  electrothennal  flow  outwards  away  from  the  constriction  opposes 
the  direction  of  bio-particle  trapping  force  under  positive  dielectrophoresis  towards  the 
constriction  tip  [9].  Hence,  we  analyze  the  influence  of  device  scaling;  i.e.  the  enhancement  of 
constriction  ratio  (R)  to  enable  trapping  of  smaller  bio-particles  and  the  reduction  of  channel 
depth  (t)  to  reduce  the  Joule  heating  on  the  net  balance  of  dielectrophoretic  trapping  force  (FDep) 
versus  electrothermal  drag  force  (Fet)  on  bio-particles.  This  balance  of  force  will  limit  the  range 
of  electric  fields  and  media  conductivity  that  can  be  applied  towards  trapping  nanoscale  bio¬ 
particles  (a  ~  5-50  mn),  where  large  constriction  ratios  are  required  as  per  the  results  presented  in 


Figure  2.  It  should  be  noted  that  the  dissipative  action  under  Fet  is  more  long-range  than  the 
trapping  action  under  FDep,  due  to  the  dependence  of  the  latter  on  field  gradients  which  are 
highly  localized  at  the  constriction  tip.  Furthermore,  while  Fdep  acts  almost  instantaneously  on 
the  polarized  particles,  Fet  requires  a  significant  temperature  gradient  to  set-in,  which  can  take  a 
longer  time  to  reach  steady  state.  These  effects  are  apparent  in  Figure  5,  which  presents 
observations  on  the  gradual  dissipation  of  dielectrophoretic  trapping  under  PDEP  due  to  the 
ensuing  electrothennal  drag  force  on  —50  mn  silica  nanoparticles  within  conductive  media 
(crm=0.5  S/m)  and  high  fields  (-350  V/cm).  The  peak  dielectrophoretic  trapping  is  reached  within 
-10  seconds  (Figure  5a),  at  which  point  the  dissipative  action  of  FEt  lowers  the  degree  of  DEP 
trapping  (Figure  5b-d).  As  per  the  fluorescence  profiles  in  Figure  5e,  beyond  the  peak  DEP 
trapping  time,  electrothermal  drag  force  causes  the  fluorescence  intensity  to  gradually  drop  and 
the  peak  intensity  is  shifted  closer  to  the  left  or  right  constriction  tips,  where  Fdep  is  highest,  and 
away  from  the  center  of  the  constriction  closer  to  the  hotspot,  where  Fet  is  highest. 

As  discussed  previously,  reducing  the  constriction  channel  depth  (t)  lowers  the  Joule 
heating  and  temperature  of  the  hotspot  at  the  center  of  the  constriction  (x,  y,  z  =  0,  0,  0)  (Figure 
4c  and  4d).  However,  in  order  for  this  to  cause  a  lowering  of  FEt,  the  ensuing  drop  in 
temperature  difference  between  the  hotspot  and  channel  walls  (which  are  close  to  ambient) 
should  be  sharper  than  the  drop  in  channel  depth  (t),  since  FEt  depends  on  the  temperature 
gradient  (VT=temperature  difference  /  depth).  It  is  apparent  from  Figure  6,  within  the  plots  of 
Fet  versus  applied  field  (Figure  6a)  and  versus  om  (Figure  6b),  that  FEt  is  lowered  with  channel 
depth,  especially  at  fields  higher  than  200  V/cm  within  media  of  am  =  1  S/m  (Figure  6a);  and  at 
am  >  1  S/m  at  E=400  V/cm  fields  (Figure  6b).  However,  it  must  be  noted  that  this  drop  in  FEt 
with  depth  is  less  sharp  than  the  drop  with  depth  of  the  respective  temperature  values  at  the 


hotspot  in  Figure  4c  and  d,  since  FEt  depends  VT.  The  greater  degree  of  lowering  of  FET  at 
successively  lower  channel  depths  (t)  is  apparent  within  the  plot  of  FEE  versus  channel  depth  in 
the  inset  of  Figure  6b,  where  FEE  exhibits  an  exponential  reduction  at  sub- urn  channel  depths, 
suggesting  that  significant  gains  to  the  balance  of  FdEp  versus  Fit  can  be  realized.  It  is  also 
apparent  that  FEE  depends  non-linearly  on  the  field  (Figure  6a)  and  linearly  on  am  (if  Figure  6b 
were  plotted  on  a  linear  X-axis),  as  expected  from  its  dependence  on  amE2. 

Finally,  we  explore  the  trade-offs  on  the  balance  of  FdEp  versus  FEE  for  the  trapping  of 
successively  smaller  bio-particles  using  successively  higher  constriction  ratios  (R).  While  the 
dielectrophoretic  forces  exhibit  a  cube-fold  dependence  on  particle  radius,  the  electrothermal 
drag  force  falls  linearly  with  particle  radius,  as  given  by:  FDrag  —  6nrjav,  where  tj  is  the  fluid 
viscosity,  a  is  the  particle  radius  and  v  is  the  fluid  velocity.  Figure  7a  summarizes  this  force 
balance  in  terms  of  the  respective  force  magnitudes  and  Figure  7b  summarizes  it  in  terms  of 
velocity  magnitudes  due  to  the  respective  forces,  to  eliminate  the  effect  of  FEj  on  particle  size. 
For  a  given  bio-particle  size,  both  the  dielectrophoretic  trapping  forces  and  the  electrothennal 
drag  forces  are  observed  to  increase  equivalently  with  increasing  constriction  ratios.  However, 
with  decreasing  bio-particle  size,  FDEp  falls  more  sharply  than  FET  due  to  its  dependence  as  cube¬ 
fold  rather  than  linear  on  particle  size,  thereby  reducing  the  gap  between  the  two  forces  at  a 
given  constriction  ratio  (R),  and  this  is  especially  apparent  at  the  smaller  bio-particle  sizes. 
Hence,  while  at  a  =  50  mn,  FEE  values  are  insignificant  due  to  the  greater  than  4  orders  of 
magnitude  lower  FET  versus  FDEp,  at  a  =  5nm,  the  difference  between  the  two  forces  is  just  two 
orders  of  magnitude  at  am=l  S/m  and  E=400  V/cm,  and  the  difference  drops  to  less  than  one 
order  of  magnitude  at  am=10  S/m  and  E=400  V/cm.  Note  that  this  analysis  is  based  on  a 
comparison  of  the  maximum  values  of  FEj  to  FdEp.  Hence,  this  difference  between  the  two  forces 


is  likely  to  be  further  lowered  if  we  had  included  the  following  within  our  analysis:  (1)  The  more 
localized  profile  of  Fdep  at  the  constriction  tip  in  comparison  to  Fet,  due  to  the  dependence  of  the 
FDep  on  field  gradient;  (2)  The  fall  of  FDep  within  media  of  higher  conductivity  (cm),  due  to 
changes  in  polarizability  with  decreasing  (ap-am)  for  the  Clausius-Mossotti  factor  (KCm  in  Eq.  1); 
(3)  The  lowered  FDep  due  to  the  limited  polarizability  of  bio-particles,  such  as  ss-DNA  and 
proteins,  to  result  in  Kcm  below  the  unity  value  assumed  within  the  calculations  herein. 

CONCLUSIONS:  In  this  work  we  investigate  the  effect  of  scaling  the  constriction  device 
geometry  on  the  balance  of  dielectrophoretic  trapping  versus  electrothennal  drag  force  on 
nanoscale  bio-particles  (<  100  mn).  Effective  dielectrophoretic  trapping  of  nanoscale  bio¬ 
particles  requires  high  constriction  ratios;  i.e.  at  least  500x  for  -50  mn  particles  and  5000x  for  ~5 
mn  particles.  However,  within  media  of  conductivities  in  the  1-10  S/m  range,  Joule  heating 
causes  temperature  rise  to  several  tens  of  degrees  above  ambient  in  the  hotspot  region  at  the 
midway  depth  point  within  the  constriction  region  and  the  resulting  electrothennal  flow  can 
cause  a  significant  dissipation  of  the  trapping  of  nanoscale  bioparticles.  Through  reduction  of  the 
depth  of  the  highly  constricted  channels  to  sub-micron  sizes,  the  degree  of  Joule  heating  can  be 
substantially  reduced  thereby  significantly  enhancing  the  range  of  applied  voltages  and  media 
conductivity  that  can  be  applied  to  enable  pre-concentration  of  silica  nanoparticles  (~50  mn)  and 
strep tavidin  protein  biomolecules  (~5  mn). 
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Figure  1(a)  Device  geometry  for  the  insulator  constriction  in  micro  fluidic  channel,  (b) 
Dimensions  and  coordinate  system  of  the  microchannel  for  the  device  modeling  studies. 
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Figure  2:  DEP  force  (a-d)  and  velocity  profiles  for  constriction  ratios  of  500x  (a-b  and  e-f, 
respectively)  and  5000x  (c-d  and  g-h)  for  trapping  particles  of  hydrodynamic  radius:  a  =  50  nm 
(a,  c,  e  and  g)  and  a  =  5  nm  (b,  d,  f,  and  h). 
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Figure  3:  (a)  Fluorescence  intensity  versus  time  for  pre-concentration  of  fluorescently  labeled 
(Alexa  488)  streptavidin  protein  molecules  in  PBS  buffer  (-150  inM  salt)  using  devices  with 
constriction  ratios  of:  (b)  9000x  and  (c)  3000x.  The  fluorescence  intensity  values  corresponding 
to  1  ug/mL  and  1  mg/mL  concentration  levels  are  shown.  The  device  shows  protein  pre¬ 
concentration  by  positive  dielectrophoresis  at  100  kHz  and  fields  of  -300  Vpp/cm.  The  channel 
depth  is  200  mn  for  both  the  constriction  ratios. 
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Figure  4:  (a)  Temperature  profile  across  channel  depth  within  the  constriction  region  shows  a 
hotspot  midway  through  the  channel,  (b)  Destruction  of  PDMS  constriction  tips  due  to  Joule 
heating  under  -400  V/cm  fields  with  am=0.5  S/m.  The  sharp  rise  in  temperature  at  the  hotspot 
with  varying:  (c)  media  conductivity  (am);  and  (d)  applied  field,  can  be  lowered  through 
reduction  in  channel  thickness  (t),  especially  within  5000x  constriction  ratio  devices. 
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Figure  5:  Dissipation  of  DEP  trapping  due  to  electrothermal  flow  is  apparent  within:  (a)  -  (d) 
fluorescence  images  of  the  constriction  tip  region  after  peak  DEP  accumulation  (~  10s);  and:  (e) 
the  fluorescence  intensity  profiles  near  the  constriction  tip. 


Electric  Field  (V/ cm) 


Figure  6:  Rise  of  electrothennal  drag  force  (FET):  (a)  quadratically  with  applied  field  (E),  and  (b) 
linearly  with  media  conductivity  (om),  can  be  reduced  through  lowering  the  channel  depth  (t). 
However,  this  reduction  is  less  sharp  than  that  of  Joule  heating  with  channel  depth  in  Figure  4c 


and  d,  due  to  the  dependence  of  FEt  on  temperature  gradient  (VT).  An  exponential  reduction  in 
Fet  is  apparent  within  the  inset  of  (b)  for  sub -pm  channel  depths. 


Figure  7:  Balance  of  dielectrophoretic  trapping  forces  versus  electrothermal  drag  forces  as  a 
function  of  particle  radius  (a),  constriction  ratio  (R)  and  media  conductivity  (am),  expressed  as: 


(a)  force  on  particles,  and  (b)  velocity  of  particles  under  this  force. 
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